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In this paper we prove that periodic boundary-value problems (EVPs) for delay 
differential equations are locally equivalent to finite-dimensional algebraic 
systems of equations. We rely only on regularity assumptions that follow those 
of the review by Hartung et al. (2006). Thus, the equivalence result can 
be applied to differential equations with state-dependent delays (SD-DDEs), 
transferring many results of bifurcation theory for periodic orbits to this class of 
systems. We demonstrate this by using the equivalence to give an elementary 
proof of the Hopf bifurcation theorem for differential equations with state- 
dependent delays. This is an alternative and extension to the original Hopf 
bifurcation theorem for SD-DDEs by Eichmann (2006). 

1. Introduction 

If a d5mamical system is described by a differential equation where the derivative at the 
current time may depend on states in the past one speaks of delay differential or, more 
generally, functional differential equations (FDEs) . A reasonably general formulation of an 
autonomous dynamical system of this type looks like this: 

x(f)=/(x„ja) (1) 

where t > is an upper bound for the delay. On the right-hand side / is a functional, 
mapping C°([— t, 0]; M") (the space of continuous functions on the interval [— t,0] with 
values in M") into E". The dependent variable x is a function on [—T:,T^ax) for some 
Tjnax > 0, and is the current /uncf ion segment: x^(s) — x{t +s) for 5 e [— t, 0] such that 
e C°([— T,0];E"). The second argument /i e W is a system parameter. For a system of 
the form (1) one would have to prescribe a continuous function x on the interval [— t,0] 
as the initial value and then extend x toward time T^^^ (see textbooks on functional 
differential equations such as [4, 11, 22]). 



A long-standing problem with certain types of FDEs is that they do not fit well into the 
general framework of smooth infinite-dimensional dynamical system theory. The problem 
occurs whenever the functional / invokes the evaluation operation in a non-trivial way, 
that is, for example, if one has a state-dependent delay. A prototypical caricature example 
would be the functional 

f:Uxm^m, /(x,jLi) = pi-x(-x(0)), (2) 

where U = {x G C°([-t,0];E) : < x(0) < t] is an open set in C°([-t,0];E). The 
corresponding FDE is 

xit)^ fi- x{t - x{t)). (3) 

Here, / evaluates its first argument x at a point that itself depends on x. We restrict 
ourselves to solutions x of (3) with x(t) e (0, t) for t > to avoid problems with causality 
and to limit the maximal delay to t (always keeping x^ in [/). 

The difficulty with (3) stems from the fact that / as a map is only as smooth as its 
argument x. Specifically, the derivative of / with respect to its first argument in this 
example exists only for x e C^([— t, 0]; M) (the space of all continuously differentiable 
functions on [— t, 0]): 

d^f : CH[-T,0];M) x M x CH[-^,0];M) ^ E, 

(A) 

d'fix,^,y)^x'i-xmy{0)-yi-xm. 

So, if we choose U as the phase space for initial-value problems (IVPs) in example (3) then 
the functional / is not differentiable for all elements oi U. In fact, it is not even locally 
Lipschitz continuous in U. Indeed, Winston [25] gave an example of an initial condition 
in U for (3) (with = and t > 1), for which the IVP did not have a unique solution. 
This counterexample is not surprising since the right-hand side / does not fit into the 
framework that the textbooks [4, 11, 22] assume to be present. A result of Walther [24] 
rescues IVPs with state-dependent delays (such as (3)) by restricting the phase space in 
general to the closed submanifold Q of C^([— t, 0]; M"): 

Q = {x e C\[-T,0];E") : x^(0) =/(x)}. 

Walther [24] could prove the existence of a semiflow inside this manifold that is continu- 
ously differentiable with respect to its initial conditions. However, this result is restricted to 
a single degree of differentiability. Results about higher degrees of smoothness are lacking 
for the semiflow [12]. 

A tjApical task one wants to perform for problems of type (1), or example (3), is bi- 
furcation analysis of equilibria and periodic orbits. Equilibria are solutions x of (1) that 
are constant in time, and periodic orbits are solutions x e C^(]R;]R") of (1) that satisfy 
x(t + r) = x(t) for some T > and all t e E. Equilibria of the general FDE (1) can be 
determined by finding the solutions (p, /i) e M" x W of the algebraic system of equations 

0=/(£oP,/i) (5) 



2 



where Eq is the trivial embedding 

£o:K"^C°([-T,0];M"), [£oP](5)=P for alls e [-t,0]. 

We observe that, even though the FDE (1) is an infinite-dimensional system, its equilibria 
can be found as roots of the finite-dimensional system (5) of algebraic equations. Moreover, 
the regularity problems of the semiflow do not affect (5): in the example (3), the algebraic 
equation (5) reads — jJL — p, which is smooth to arbitrary degree, and can be solved even 
for negative jJL (near equilibria with /i = p < the semiflow does not exist). 

In this paper we establish a system similar to (5), but for periodic orbits: we find a 
finite-dimensional algebraic system of equations that does not suffer from the regularity 
problems affecting the semiflow, and an equivalence between solutions of this algebraic 
system and periodic orbits of (1). In comparison, for ordinary differential equations (ODEs) 
of the form x(t) = /(x(t), pt) with a smooth / : E" x E", the fact that the problem 
of flnding periodic orbits can be reduced to algebraic root-flnding is well known [9]. For 
example, in ODEs one can use the algebraic system = X{T; p,ix)—p where t X{t; p, ji) 
is the trajectory deflned by the IVP starting from p e E" and using parameter /i e . 

A central notion in the construction of the equivalent algebraic system for periodic orbits 
of FDEs are periodic boundary-value problems (EVPs) for FDEs on the interval [—n, n] with 
periodic boundary conditions (which we identify with the unit circle T) . Periodic orbits of 
(1) can then be found as solutions of periodic EVPs. If one wants to make the equivalence 
result useful in practical applications, one has to find a regularity (smoothness) condition on 
the right-hand side / that includes the class of state-dependent delay equations reviewed 
in [12], while still ensuring that it is possible to prove the existence of an equivalent 
algebraic system. We use exactly the same condition as used by Walther in [24] to prove 
the existence of a continuously differentiable semiflow, the so-called extendable continuous 
differentiability (originally introduced as "almost Frechet differentiability" in [21]), which 
implies a restricted form of local Lipschitz continuity. We generalize restricted continuous 
differentiability to higher degrees of restricted smoothness (which we call EC^ smoothness) 
in a similar fashion as Krisztin [19] did for the proof of the existence and smoothness of 
local unstable manifolds of equilibria. Our deflnition of EC'^ smoothness is comparatively 
simple to state and check, and lends itself easily to inductive proofs. 

After introducing the notation for periodic EVPs and EC'^ smoothness we state the 
main result, an equivalence theorem between periodic EVPs and algebraic systems of 
equations in Section 2. The equivalence theorem reduces statements about existence and 
smooth dependence of periodic orbits of FDEs to root-finding problems of smooth algebraic 
equations. The result is weaker than the corresponding results for equilibria of FDEs and for 
periodic orbits of ODEs because the equivalence is only valid locally. For any given periodic 
function Xg with Lipschitz continuous time derivative we construct an algebraic system that 
is equivalent to the periodic EVP in a sufficiently small open neighborhood of Xq. However, 
the result is still useful, as we then demonstrate in Section 3. We apply the equivalence 
theorem in the vicinity of equilibria for which the linearization of (1) has eigenvalues on 
the imaginary axis (for example, near Xq = jj, = 7i/2 in example (3)) to prove the Hopf 
Eifurcation Theorem. The equivalence theorem reduces the proof of the Hopf Eifurcation 
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Theorem to an application of the Algebraic Branching Lemma [1]. This provides a complete 
proof for the Hopf Bifurcation Theorem for FDEs with state-dependent delays, including the 
regularity of the emerging periodic orbits. We discuss differences to the first version of the 
proof by Eichmann [5] and the approach of Hu and Wu [13] in Section 3. The equivalence 
is applicable in other scenarios where one would expect branching of periodic solutions. 
Examples are period doublings, the branching from periodic orbits with resonant Floquet 
multipliers on the unit circle in Arnol'd tongues, and branching scenarios in FDEs with 
symmetries. We give a tentative list of straightforward applications and generalizations of 
the equivalence theorem in the conclusion (Section 4). 

We note that the theorem stated in Section 2 differs from statements about numerical 
approximations. As part of the theorem we also provide a map X that maps the root of 
the algebraic system back into a function space to give the exact solution of the periodic 
BVI^ and a projection P that maps functions to finite-dimensional vectors (and, hence, 
periodic orbits to roots of the algebraic system). In numerical methods one typically 
has to increase the dimension of the algebraic system in order to get more and more 
accurate approximations of the true solution whereas the dimension of the algebraic system 
constructed in Section 2 is finite. 



This section states the assumptions and conclusions of the main result of the paper, the 
Equivalence Theorem stated in Theorem 2.5. Before doing so, we introduce some basic 
notation (function spaces on intervals with periodic boundary conditions and projections 
onto the leading Fourier modes) . 

Periodic EVPs 

We first state precisely what we mean by periodic BVP and introduce the usual hierarchy 
of continuous, continuously differentiable and Lipschitz continuous functions on the 
compact interval [ — 7i, 7i] with periodic boundary conditions. For j > we will use the 
notation C-'(T; E") for the spaces of all functions x on the interval [—n, n] with continuous 
derivatives up to order j (including order and j) satisfying the periodic boundary 
conditions x'^'^(— 7i) = x'^'^\n) for I = 0...j. Elements of C°(T;E") are continuous and 
satisfy x(— 7i) = x(7i). For derivatives of order j > 0, x^^\—7i) is the right-sided jth 
derivative of x in —n, and x'^^\7i) is the left-sided jth derivative of x in n. The norm in 
C^'(T;E") is 



We can extend any function x in C-'(T; M") to arguments in R by defining x(t) = x(t — 2kn) 
where k is an integer chosen such that —n < t — 2kn < n (we will write t^^^^-n.-a) later). 
Thus, every element of C^ (T; E") is also an element of BC-' (E; E"), the space of functions 
with bounded continuous derivatives up to order j on the real line. We use the notation 
t e T for arguments t of x, and also call T the unit circle. This make sense because the 



2. The Equivalence Theorem 
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parametrization of the unit circle by angle provides a cover, identifying T with E where we 
use [—71, 7i) as the fundamental interval. 

Additional useful function spaces are the space of Lipschitz continuous functions and, 
correspondingly, spaces with Lipschitz continuous derivatives, denoted by C-''^(T;E"), 
which are equipped with the norm 

J |xW(5)-xW(f)l 

\\x\\jj — max < \\x\\j, sup ■ ■ 

I t,seK 1^ ~ t| 

(x'^°^(t) refers to x(f)). Note that we used the notation t,s e M in the index of the 
supremum, as we can apply arbitrary arguments in E to a function x e C°(T;M") by 
considering it as an element of BC°(E; M"), as explained above. We use the same notation 
(C-'(J;E") and C-''HJ;E")) also for functions on an arbitrary compact interval J c E 
without periodic boundary conditions (and one-sided derivatives at the boundaries). As 
any function x e C^(T;E") is also an element of BC-'(M;M"), it is also an element of 
C-'(J;E") for any compact interval J (and the norm of the embedding operator equals 
unity). On the function spaces C-'(T;]R") we define the time shift operator 

A, : C^(T;E") ^ C-'(T;E"), [A^xKs) = x(t +s). 

The operator A^. is linear and has norm 1 in all spaces C^(T;M"). Similarly, A;, maps also 
C^'HT; M") ^ C-'''^(T; E"), and has unit norm there as well. 
Let / be a continuous functional on the space of continuous periodic functions, that is, 

/ : C°(T;E")^M". 

The right-hand side /, together with the shift A^, creates an operator in C°(T; E"), defined 
as 

F:C°(T;E")^C°(T;E") [F(x)](t) = /(A,x). (7) 

The operator F is invariant with respect to time shift by construction: F(AfX) = AfF(x). 
We consider autonomous periodic boundary-value problems for differential equations 
where / is the right-hand side: 

x(t)=/(A,x) = F(x)(t). (8) 

A function x e C^(T;E") is a solution of (8) if x satisfies equation (8) for all f e T (for 
each t e T equation (8) is an equation in E"). In contrast to the introduction we do not 
expressly include a parameter pi as an argument of / . This does not reduce generality 
as we will explain in Section 3. The main result, the Equivalence Theorem 2.5, will be 
concerned with equivalence of the periodic BVP (8) to an algebraic system of equations. 
The notion of the shift A,, on the unit circle and the operator F, combining / with the shift, 
is specific to periodic EVPs such that the BVP (8) looks different from the IVP (1) in the 
introduction. Several results stating how regularity of / transfers to regularity of F are 
collected in Appendix A. 
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Definition of EC'^ smoothiness and local (restricted) EC Lipschitz continuity 

Continuity of the functional / is not strong enough as a condition to prove the Equivalence 
Theorem. Rather, we need a notion of smoothness for / . However, as explained in the 
introduction, we cannot assume that / is continuously differentiable with degree > 1, if 
we want to include examples such as f{x) — — x(— x(0)) (see FDE (3) for = 0) into the 
class under consideration. 

The review by Hartung et al. [12] observed the following typical property of functionals 
/ appearing in equations of type (8): the derivative d^f{x) of / in x as a linear map from 
C^T; E") into E" can be extended to a bounded linear map from C°(T; E") into E", and 
the mapping 

: CHT;E") X C°(T;E")^E" defined by (x,y) ^ 5 V(x,y) 

is continuous as a function of both arguments. In other words, the derivative of / may 
depend on x' but not on y'. For the example /(x) = — x(— x(0)) this is true (see (4)). 
Most of the fundamental results establishing basic dynamical systems properties for FDEs 
with state-dependent delay in [12] rest on this extendability of d^f . 

We also rely strongly on this notion of extendable continuous differentiability. The precise 
definition is given below in Definition 2.1. In this definition we permit the argument range 
J to be any compact interval or T. We use the notation of a subspace of higher-order 
continuous differentiability not only for C^ (J;E") but also for products of such spaces in a 
natural way. Say, if 

D^C'''iJ;R"'')x ...xC^'{J;R"'^), (9) 

where I > 1, and kj > and nij > 1 are integers, and denoting the natural maximum norm 
on the product D by 

\\x\\D = \\ixi,...,Xf)\\o^ max ||x.||.., 

je{l,...,t} ' 

then for integers r > the space is defined in the natural way as 

= C'^i+''(J;E'"i) X ... X C^t+'-iJ-R"'^), with 
||x|L r = max HxPlL.. 

' 0<j<r ' 
l<i<t 

For the simplest example, D = C°(J; E"), is C^{J; E"). If J = T then the time shift 
extends naturally to products of spaces: 

AfX = (AfX^, . . . , AfXf ) for x = (x^, . . . , x^) e D. 

Definition 2.1 (Extendable continuous differentiability EC^) 

Let D be a product space of the type (9), and let / : D E" be continuous. We 
say that / has an extendable continuous derivative if there exists a map d 

I d^f -.D^xD^R" 
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that is continuous in both arguments (u, v) e x D and hnear in its second 
argument v e D, such that for all u e 




\f(u + v)-f(u)-djiu,v)\ 

Md,1 



= 0. 



(10) 



We say that f is k times continuously differentiable in this extendable sense if 
the map d^f, recursively defined as d^f — d^[d^~^f], exists and satisfies the 
limit condition (10) for d^~^f. We abbreviate this notion by saying that / is EC^ 
" smooth in D. 

The limit in (10) is a limit in E. For k — 1 the definition is identical to property (S) in 
the review [12], one of the central assumptions for fundamental results on the semiflow. 
Extendable continuous differentiability requires the derivative to exist only in points in 
and with respect to deviations in D^, but it demands that the derivative must extend in its 
second argument to D {d is linear in its second argument) . This is the motivation for 
calling this property extendable continuous differentiability. 

The definition of EC'^ smoothness for fc > 1 uses the notation that a functional (say, d^f) 
of two arguments (say, uG and v e D) for which one would write d^f{u, v), is also a 
functional of a single argument w = (u, v) e x D, such that one can also write d^f{w). 
When using this notation we observe that the space x D is again a product of type (9) 
such that d is again a functional of the same structure as / . For example, let us consider 
the functional / : x — > — x(— x(0)) from example (2) (setting fi — 0). The functional is 
well defined and continuous also on D = C°(T;E). Moreover, / is EC*^ smooth in D to 
arbitrary degree k. Its first two derivatives are: 



As one can see, the first derivative d has the same structure as / itself if we replace 
D - C°(T;E") by x D. So, it is natural to apply the definition again to 5 ^ on the space 
X D. 

Assuming that / is EC^ smooth on C°(J; E") implies classical continuous differentiability 
of / as a map from C^(J; E") into E" and is, thus, strictly stronger than assuming that / is 
continuously differentiable on CnTr;E"). 

Since every element of C^(T; E") is also an element of C^{J; E") for any compact interval 
J (and the embedding operator has unit norm), any EC^ smooth functional / : C°(J; E'") — > 
E" is also a EC'' smooth functional from C°(T; E'") into E". 

It is worth comparing Definition 2.1 with the definition for higher degree of regularity 
used by Krisztin in [19]. With Definition 2.1 the kth derivative has 2^ arguments. In 
contrast to this, the fcth derivative as defined in [19] has only k + 1 arguments (the 



: CHT;E) X C°(T;E), 
d^fiu,v) = u^(-u(0))v(0) - v(-u(0)), and 
d^f : [C^(T;E) x C\T;E)] x [c\T;E) x C°(T;E)] , 
82/ iu, V, w, x) = - u '(-u(0)) w(0) v(0) + u (-u(O)) x(0) 



(11) 



(12) 



+ w'i-uiO)) v(0) - v'(-u(0)) w(0) - x(-u(0)). 
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first argument is the base point, and the derivative is a fc-linear form in the other k 
arguments). The origin of this difference can be understood by looking at the example 
fix) — — x(— x(0)) and its derivatives in (11)-(12). Krisztin's definition applied to the 
second derivative does not include the derivative of d with respect to the linear second 
argument v (as is often convention, because it is the identity). One would obtain the 
second derivative according to Krisztin's definition by setting x = in (12). Indeed, 
the terms containing the argument x in (12) are simply d^f{u,x), as one expects when 
differentiating d^f(u,v) with respect to v, calling the deviation x. While in practical 
examples it is often more economical to use the compact notation with fc-forms, inductive 
proofs of higher-order differentiability using the full derivative only require the notion of 
at most bi-linear forms, making them less complex. 

If / is EC^ smooth then it automatically satisfies a restricted form of local Lipschitz 
continuity [12], which we call local EC Lipschitz continuity: 

Definition 2.2 (Local EC Lipschitz continuity) 

I We say that / : C°(T; E") R" is locally EC Lipschitz continuous if for every 
Xq e C^(T;E") there exists a neighborhood [/(xq) c CHT;E") and a constant K 
- such that 

I \fiy)-fiz)\<K\\y-z\\, (13) 

holds for all y and z in U{xq). 

That EC^ smoothness implies local EC Lipschitz continuity has been shown, for example, 
in [24] (but see also Lemma A.3 in Appendix A). Note that the estimate (13) uses the 
II • llo-norm for the upper bound. This is a sharper estimate than one would obtain using 
the expected || • ||i-norm. The constant K may depend on the derivatives of the elements in 
L/(xo) though. For example, for /(x) = — x(— x(0)) as in (2) with = 0, one would have 
the estimate 

+ < [l + l|x'|lo] llyllo suchthat iC<l+ max llxlli. 

This means that in this example, the neighborhood i7(xo) can be chosen arbitrarily large as 
long as it is bounded in C^(T; E"). 

The following lemma states that we can extend the neighborhood U{x) in Definition 2.2 
into the space of Lipschitz continuous functions (C°'^ instead of C^) and include time shifts 
(which possibly increases the bound K) . 

Lemma 2.3 {EC Lipschitz continuity uniform in time) 

Let / be locally EC Lipschitz continuous, and let Xg be in C°'^(T; E"). Then there 
exists a bounded neighborhood U{xq) c C°'^(T; M") and a constant K such that 

I |/(A,y) -/(A,z)| < K\\A,y - A,z||o = K\\y - z\\o 

holds for all y and z in f/(xo), and for all t e T. Thus, ||F(y) - F(z)||o < if || j -z||o 
for all y and z in U{xq). 
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Recall that F(x)(t) — f(A^x). See Lemma A.3 and Lemma A.5 in Appendix A for the proof 
of Lemma 2.3. A consequence of Lemma 2.3 is that the time derivative of a solution Xq 
of the periodic BVP is also Lipschitz continuous (in time): if Xo(f ) = / (A^Xq) then there 
exists a constant K such that 



\x'^it)-x'^is)\\,<K\t-s\ 



(14) 



Thus, Xq e C^'^(T;E"). This follows from Lemma 2.3 by inserting A^Xq and A^Xg for y 
and z and using that x'^it) = /(A^Xq) (it is enough to show (14) for \t —s\ small). 

Projections onto subspaces spanned by Fourier modes 

The variables of the algebraic system in the Equivalence Theorem will be the coefficients 
of the first N Fourier modes (where N will be determined as sufficiently large later) of 
elements of C°'^(T;E") (the space of Lipschitz continuous functions on T). Consider the 
functions on T 



bo — t ^ -, bj^ — t ^ cos(fct), b_j^ = t >-> sin(fct) 

for = 1, . . . , GO (which is the classical Fourier basis of L^(T; M)). For any m > 1 we define 
the projectors and maps 



P^:C^(T;E'")-C^(T;E'"), [P^x](f); = ^ 



1 

71 



Wit), 



I-Pn, 



N 



: C^(T;E'") ^ E-x(2w+i)^ [-RjvxJi^fc 
L : C^(T;E'")^ C^(T;E'"), [Lx](t) 



k=-N 

n 



(15) 



fefc(5)Xi(5)ds, 



x(5) —RqX ds — 



Qo[x](5)d5. 



The projector Pj^ projects a periodic function onto the subspace spanned by the first 2N + 1 
Fourier modes, and Qjy is its complement. The map £^ maps a vector p of 2N + 1 Fourier 
coefficients (which are each vectors of length n themselves) to the periodic function that 
has these Fourier coefficients. The map R^^ extracts the first 2N + 1 Fourier coefficients 
from a function. The simple relation P^^ — Ej^Rj^ holds. The vector RqX is the average of a 
function x, and Qg subtracts the average from a periodic function. The operator L takes 
the anti-derivative of a periodic function after subtracting its average (to ensure that L 
maps back into the space of periodic functions) . In all of the definitions the degree j of 
smoothness of the vector space can be any non-negative integer. The operator L not 
only maps back into itself, but it maps C-''(T; R"") into C-'+HT; M'"). 
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We do not attach an index m to the various maps to indicate how many dimensions 
the argument and, hence, the value has because there is no room for confusion: for 
example, if x e C°(T; E^) then Pf^x e C°(T; E^) such that we use the same notation Pj^x 
for X : T '-^ E"" with arbitrary m. Similarly, we apply all maps also on product spaces D of 
the type C'^^iT; x . . . x C^^'CT; ^""0 introduced in Equation (9) by applying the maps 
element-wise. For example, 

Pf^x = (Pjv^i, ■ ■ ■ , Pn^i) for X = (xi, ...,Xi)GD, 

E^p = (£^pi, . . .,£jvP£) for p = (Pi, . ..,Pt) e E'"i^(2iv+i) ^ ^ ^m,x(2jv+i)_ 



Equivalent integral equation 

We note the fact that a function x e C^T; E") solves the periodic BVP x(t) = /(A^x) = 
F(x)(t) if and only if it satisfies the equivalent integral equation 



x(t) = ^(0) + 



F(x)(s)ds for all t ST. (16) 



JO 



For each t e T, Equation (16) is an equation in E". In particular, the term x(0) is in E". 
Thus, the integral equation (16) is very similar to the corresponding integral equation 
used in the proof of the Picard-Lindelof Theorem for ODEs [3] . This is in contrast to the 
abstract integral equations used by Diekmann et al. [4] to construct unique solutions to 
IVPs, in which equality at every point in time is an equality in function spaces. It is the 
similarity of (16) to its ODE equivalent that makes the reduction of periodic EVPs to finite 
dimensional algebraic equations possible. One minor problem is that the Picard iteration 
for (16) cannot be expected to converge. In fact, the integral term F(x)(s)ds does not 
even map back into the space C°(T; E") of periodic functions, even if x is in C°(T; E"). 
However, a simple algebraic manipulation using the newly introduced maps L, P^, Q^, 
and Rj^ removes this problem (remember that F(x)(t) = /(A^x)): 

Lemma 2.4 (Splitting of BVP) 

I Let N > be an arbitrary integer. A function x e C°(T; E") and a vector p e 

^nx(2]V+l) g^^jgfy 

x(t)=/(AfX) and p^R^x, (17) 

if and only if they satisfy the system 

x = FwP+QjvLF(x), (18) 
= [PoFix) + Qo [E^p - P^LFix)) ] . (19) 



Note that the map Rf^ extracts the lowest 2Ar-l-l Fourier coefficients from a periodic function. 
Equation (18) can be viewed as a fixed-point equation in C°'^(T; E"), parametrized by p. 
We will apply the Picard iteration to this fixed-point equation instead of (16). Equation (19) 
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is an equation in ]R"><C2w+ij the Picard iteration converges then the fixed-point equation 
(18) can be used to construct (for sufficiently large N) a map X : U <z ]R"><(2]v+i) ^ 
C°'^(T;E"), which maps the parameter p to its corresponding fixed point x. Inserting this 
fixed point x —X{p) into (19) turns (19) into a system of n x {2N + 1) algebraic equations 
for the n X {2N + l)-dimensional variable p, making the periodic BVP for x equivalent to 
an algebraic system for its first 2N + 1 Fourier coefficients, p. The proof of Lemma 2.4 is 
simple algebra, see Section 5.2. 

Statement of the Equivalence Theorem 

Using the Splitting Lemma 2.4 we can now state the central result of the paper. The 
intention to treat (18) as a fixed-point equation motivates the introduction of the map 

Mjv : C°'\T; M") x M"x(2w+i) ^ c°'HT; M") given by 
M^{x,p)^E^p + Q^LFix). 

This means that we will look for fixed points of the map Mjv(-,p) for given p and sufficiently 
large N. We will do this in small closed balls in C°'^(T;]R") (the space of Lipschitz 
continuous functions) such that it is useful to introduce the notation 

Bfix,) = {x e C°'HT;E") :\\x-x,\l^^<5], 

for 5 > and Xq e C°'HT;M"). That is, ^"'^(xo) is the closed ball of radius 5 around 
Xq e C^'HTTjK") in the || • Ho^^-norm (the Lipschitz norm on T). 

Theorem 2.5 (Equivalence between periodic BVPs and algebraic systems of equations) 

Let / be £C^™^ smooth, and let Xq have a Lipschitz continuous derivative, that is, 
Xg e C^'^(T; E"). Then there exist a 5 > and a positive integer N such that the 
map Mjv(-,p) has a unique fixed point in 5°^ (xg) for all p in the neighborhood U 
of Rj^Xq given by 

I U={pG M"x^^^+^) : ||£^ [p - R,x,] < 25} . 

The maps 

X : [/^ C°(T;E"), X(p) = fixed point of Mjv(-,p) inB°g^(xo), 
g:U^ Enx(2iv+i)^ g(p) [PoF(X(p)) + Qo [E^p - P^LFiXip)))] , 

are j^^^ times continuously differentiable with respect to their argument p, and 
Xip) is an element of C-'">^'<+^(T; E"). Moreover, for all x e B°'^(xo) the following 
equivalence holds: 

I x(f)=/(A,x) 
if and only if p = Rj^x is in U and satisfies 
I g(p) = and x=X(p). 
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Theorem 2.5 is the central result of the paper. It implies that, for any Xq e C^'^(T; E") all 
solutions of the periodic BVP in a sufficiently small neighborhood of Xq lie in the graph 
of a finite-dimensional manifold. Moreover, these solutions can be determined by finding 
the roots of g in [/ c ]r"x(2Jv+i) ^^^^ ^^lat Theorem 2.5 is different from statements 
about numerical approximations. Even though the integer N is finite, solving the algebraic 
system g(p) = and then mapping the solutions with the map X into the function space 
C°(T;E") gives an exact solution x —X{p) of the periodic BVP x(t) =/(AfX). 

The size of the radius 5 of the ball in which the equivalence holds depends on how large 
one can choose 5 such that a local EC Lipschitz constant K for F exists for 5°^ (xg) (such 
neighborhoods exist according to Lemma 2.3). In many applications (in particular, in the 
example (3)) this can be any closed ball in which the right-hand side / is well defined 
(at the expense of increasing K for larger balls). Once the local EC Lipschitz constant 
K is determined, one can find a uniform upper bound R for the norm ||F(x)||o i for all 
X e Bgg (xq) (see Lemma A.5). The integer N, which determines the dimension of the 
algebraic system, is then chosen depending on R, K and ||Xq||o j. 

Section 5 contains the complete proof of Theorem 2.5. The first step of the proof of 
Equivalence Theorem 2.5 is the existence of the fixed point of Mjv in 5°^ (xq) for p gU. 
This is achieved by applying Banach's contraction mapping principle to the map Mf^(-,p) 
in the closed ball (xg). The only peculiarity in this step is that we apply the principle 
to Bg^^(xo), which is a closed bounded set of Lipschitz continuous functions, using the 
(weaker) maximum norm (|| • ||o). This is possible because closed balls in C°'^(T; E") are 
complete also with respect to the norm || • ||o- With respect to the maximum norm the map 
Mjy(-,p) : X Ej^p + QjvLF(x) becomes a contraction for sufficiently large N (because the 
norm of the operator Qj^L is bounded by C\og(N)/N , and F has a Lipschitz constant K 
with respect to || • ||o in 5°^ (xq)). 

After the existence of the fixed point of Mjv(-,p) is established in Section 5.3 the equiva- 
lence between the algebraic system g(p) = and the periodic BVP x(t) = /(A(.x) in the 
smaller ball B°'^(xo) follows from the Splitting Lemma 2.4. 

The smoothness (in the classical sense) of the maps X and g follows, colloquially 
speaking, from implicit differentiation of the fixed-point problem x = E^^p + QjvLF(x) 
with respect to p. Section 5.5 checks the uniform convergence of the difference quotient 
in detail. Section 5.6 uses the higher degrees of EC^"^ smoothness of / to prove higher 
degrees of smoothness forX and g. For proving higher-order smoothness one has to check 
only if the spectral radius of a linear operator is less than unity, but the inductive argument 
requires more elaborate notation than the first-order continuous differentiability. 

3. Application to periodic orbits of autonomous FDEs — Hopf 
Bifurcation Tlieorem 

Let us come back to the original problem, the parameter-dependent FDE (1) x =/(Xf,/i), 
where jJL e W is a system parameter and the functional / : C°(J;E") x W E" is defined 
for first arguments that exist on an arbitrary compact interval J. Periodic orbits are solutions 
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X of i: — f(Xf, jLt) that are defined on M and satisfy x(t) — x(t + T) for some T > and all 
t e E. 

Let X be a periodic function of period T — 27i/co. Then the function y(s) = x(s/<d) is a 
function of period 27i (s e T). This makes it useful to define the map 

S :BC°(E;E") x E^BC°(E;M") [S(y, fo)](s) = y(cos), 

such that S(y,co)(f) = x(t) for all t e E (remember that BC°(M;M") is the space of 
bounded continuous functions on the real line). Then x e C^{R; E") satisfies the differential 
equation 

x(f)=/(x„pi) (20) 

on the real line and has period 27i/£o if and only if y = S(x, 1/co) e C^(T;E") satisfies the 
differential equation 

1 

y(s) = — /(S(A,y,co),pi). 

CO 

Let us define an extended differential equation 



^ext(5)=/ext(A,Xe^J, (21) 

where f^^, maps C°(T; E"+^+' ) into and is defined by 



fe: 



CO 



CUt( CO ) 



/ {Siy,RoCo),Ro^) / cut{RoCo) 





where 



CO if CO > CO cutoff > 

smooth, uniformly non-negative extension for co < co (.^off 



for y G C°(T;E"), co e C°(T;E) and pi e C°(T;M'') (recall that takes the average 
of a function on T) . We have used in our definition that any functional / defined for 
X e C°(J; E") is also a functional on C°(T; M") (periodic functions have a natural extension 
x(t) = ■x(tmod[-7r,7t)) if f s M is arbitrary). The extended system has introduced the 
unknown co and the system parameter [jl as functions of time, and the additional differential 
equations co — 0, fi — 0, which force the new functions to be constant for solutions of (21). 
We have also introduced a cut-off for co close to zero to keep /g^t globally defined. The 
extended BVP (21) is in the form of periodic EVPs covered by the Equivalence Theorem 2.5. 
Thus, if /gxt is EC^"'"' smooth then BVP (21) satisfies the assumptions of Theorem 2.5 in 
the vicinity of every periodic function Xg ext ^ C^'^(T; M""*"^"*"^). Any solution x^^t = (y, oj, />i) 
that we find for (21) corresponds to a periodic solution t yicot) of period 2ti/RqC0 at 
parameter RqIjl for (20) and vice versa, as long as RqCO > co cutoff- The condition of EC^^'^ 
smoothness has to be checked only for the first n components of the function /gxt since its 
final V -I- 1 components are zero. 

Application of the Equivalence Theorem 2.5 results in a system of algebraic equations 
that has (n-l- v -I- l)(2Ar -1-1) variables and equations, where N is the positive integer proven 
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to exist in Theorem 2.5. Let us denote as F — {Fy,F^,F^) the components of the right- 
hand side Fgxt (given by Fext(-^ext)(f ) — /extC'^f^ext))^ of which F^ and F^ are identically 
zero. Let p = iPy,Pco'Pix) be the 2N + 1 leading Fourier coefficients of y, co and jji, 
respectively (these are the variables of the algebraic system constructed via Theorem 2.5), 
and X(p) = (X^(p),X„(p),X^(p)) be the map from j?("+^'+i)(2iv+i) i^to C^--(T; 
Then several of the components of p can be eliminated as variables, and the equations for 
p resulting from Theorem 2.5 correspondingly simplified. Since F is identically zero in its 
last V + 1 components we have 

Hence, the right-hand side 

g(p) = [PoFiXip)) + Qo (£jvP - PMLFiXip))) ] , 

defined in Theorem 2.5, has v + 1 components that are identical to zero (since PqF^X^p)) = 
for the equations co = and /i = 0). Furthermore, g(p) = contains the equations 
RjvQo^nPw = ^nd R^QqEi^Pij^ = 0, which require that all Fourier coefficients (except the 
averages RqCO and RqIJl) of jJL and co are equal to zero. This means (unsurprisingly) that 
the algebraic system forces co and ji to be constant. Thus, we can eliminate Rj^QqE^^p^ 
and RjvQgEjvP/i (which are 2N{v + 1) variables), replacing them by zero, and drop the 
corresponding equations. Since co and /i must be constant, we can replace the arguments 
p^ and p^ of X by the scalar RqE^^p^ (which we re-name back to co) and the vector 
Rq^jvP/x ^ (which we re-name back to pi). This leaves the first n{2N + 1) algebraic 
equations 

=Rjv {PoFyiXyiPy, CO, ix), CO, /i) + Qo [Ej^Py - Pt^LFy{Xy{py, co, pi), co, pi)] ) , (22) 

which depend smoothly (with degree Jmax) the n{2N + 1) variables Py and the parameters 
to e E and pi e E''. Overall, (22) is a system of n x (2N -I- 1) equations. 



Rotational Invariance 

The original nonlinearity F, defined by [F(x)](t) = /(A^x) is equivariant with respect to 
time shift: AfF(x) = F(AtX) for all f e T and x e C°(T;E"). Furthermore, A^ commutes 
with the following operations: 

AfQjvL = QiviA, (if N > 0) and A^P^ = P^A^. 

This property gets passed on to the algebraic equation in the following sense: let us define 
the operation A^. for a vector p in which we consider as a vector of Fourier 

coefficients of the function £^p e C°(T; E"), by 

AfP =i?wAt£wP- 

With this definition Af commutes with Rj^ and Fjv- It is a group of rotation matrices: Af 
is regular for all t, and A2fcjr is the identity for all integers k. The definition of X(p) as a 
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fixed point of x '-^ Ej^p + Qf^LF(x) implies that AfX(p) = X(Afp). From this it follows 
that the algebraic system of equations is also equivariant with respect to . If we denote 
the right-hand-side of the overall system (22) by G{py, <D,/i) then G satisfies 

AfiiPy, 00,11) = G(AtP3„co,pi) for all t e T, p^, e ^"(^2^+1^ co > and e M'' . 
Application to Hopf bifurcation 

One useful aspect of the Equivalence Theorem is that it provides an alternative approach to 
proving the Hopf Bifurcation Theorem for equations with state-dependent delays. The first 
proof that the Hopf bifurcation occurs as expected is due to Eichmann [5]. The reduction 
of periodic boundary-value problems to smooth algebraic equations reduces the Hopf 
bifurcation problem to an equivariant algebraic pitchfork bifurcation. 
Let us consider the equation 

0=/(£oXo,Pt) (23) 

where / : C°(J; E") x E E", /i e E, J is a compact interval, Xq e E", and the operator 
£0 (as defined in (15) in Section 2) extends a constant to a function on T (and thus, on J). 
This means that (23) is a system of n algebraic equations for the n + 1 variables (xq, jU,). The 
definition of EC*^ smoothness does not cover functionals that depend on parameters. We 
avoid the introduction of a separate definition of EC*^ smoothness for parameter-dependent 
functionals that distinguishes between parameters and functional arguments. We rather 
extend Definition 2.1: 

Definition 3.1 (£C^ smoothiness for parameter-dependent functionais) 

Let J = [a, b] be a compact interval (or J = T), and D be a product space of the 
form D = C^'iJ; E'"i) x . . . x C'^'iJ; E'"0 where £ > 1, and kj > and ruj > 1 are 
integers. We say that / : D x E"^' E" is EC^ smooth if the functional 

I (x,y)eDxC°(J;E^)^/(x,j(a))eE" (24) 

is EC^ smooth (if J = T we use a = — 71). 

Requiring EC*^ -smoothness of the parameter-dependent functional / in this sense, implies 
that the algebraic system = /(Eq-^ojM) is k times continuously differentiable. Let us 
assume that the algebraic system = /(Eq-^o^M) has a regular solution Xq^jl) e E" for 
jjL close to 0. Without loss of generality we can assume that Xq^jl) = 0, otherwise, we 
introduce the new variable x„g^^ = x^i^ — EoXo(jLt) e C°(J; E"). Hence, /(O, pi) = for all jj, 
close to 0. 

The EC^ derivative of / in (0,//) is a linear functional, mapping C°(J;E"+'' ) into E". Let 
us denote its first n components (the derivative with respect to the first argument x of /) 
by a{iJL). The linear operator a(pi) can easily be complexified by defining a(/i)[x + iy] — 
a(/>i)[x] +ia(/x)[y] for x -h iy e C°([-t, 0]; C"). If / is EC*" smooth with > 2 then the 
n X n-matrix K{X,ii) (called the characteristic matrix), defined by 

if(A,p) V = Av - a(pi)[vexp(Af)] (25) 
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is analytic in its complex argument A and k — 1 times differentiable in its real argument /i 
(since the functions t ^ vexp(At) to which a(/i) is applied are all elements of C*^(J;E")). 
The Hopf Bifurcation Theorem states the following: 

Theorem 3.2 (Hopf bifurcation) 

Assume that / is EC'^ smooth {k > 2) in the sense of Definition 3.1 and that the 
characteristic matrix K(A, pi) satisfies the following conditions: 

1. (Imaginary eigenvalue) there exists an cOq > such that det K{icoQ, 0) = 
and icoQ is an isolated root of A detK{X, 0). We denote the corresponding 
null vector by = + iv; e C" (scaling it such that Iv^l^ + |v;p = 1). 

2. (Non-resonance) detK{ikco, 0) 7^ for all integers k ^ ±1. 

3. (Transversal crossing) The local root curve jJL A(pi) of detK{X,iJi) that 
corresponds to the isolated root icog at pi = (that is, A(0) = icoq) has a 
non-vanishing derivative of its real part: 

d 

07^— ReA(/i) 

pi = 

Then there exists a — 1 times differentiable curve 

p G {-£,£) ^ (x,co,pi) e CHTr;E") x E x E 
such that for sufficiently small e > the following holds: 

1. x(co-) (or S{x, co)) is a periodic orbit of x{t) — f{x^,jjL) of period 2n/co, that 
is, X e CHT;E") and 

^(f) = -/(S(A,x,co),pi), (26) 

CO 

2. the first Fourier coefficients of x are equal to (0, ^), that is. 



1 



= 

71 

1 

/3 = — 

71 



Re [viexp(it)]^x(f)df, and 



n 

en 

Im [v^ exp(it)] ^ x(t)dt, 



(27) 



3. ^1/3=0 — 0, iLi|^=o — ^i^d co\p^Q — coq, that is, the solution x, the system 
parameter jJL and the frequency co of x, which are differentiable functions of 
the amplitude 13, are equal to x — 0, fi — 0, co — cOq for 13 —0. 

The statement is identical to the classical Hopf Bifurcation Theorem for ODEs in its assump- 
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tions and conclusions apart from the regularity assumption on / specific to FDEs. Note 
that the existence of the one-parameter family (parametrized in ^) automatically implies 
the existence of a two-parameter family for ^ / due to the rotational invariance: if x is a 
solution of (26) then A^x is also a solution of (26) for every fixed s e T. Condition (27) 
fixes the time shift s of x such that x is orthogonal to Re[vi exp(it)] using the scalar 
product on T. 

The proof of the Hopf Bifurcation Theorem is a simple fact-checking exercise. We have 
to translate the assumptions on the derivative of / : C°(J; M") ^ E" into properties of the 
right-hand side of the nonlinear algebraic system (22) near (x, co, jj.) — (0, cOq, 0), and then 
apply algebraic bifurcation theory to the algebraic system. The only element of the proof 
that is specific to functional differential equations comes in at the linear level: the fact that 
the eigenvalue ico^ is simple implies that the right nuUvector e C" and any non-trivial 
left nuUvector satisfy 



This is the generalization of the orthogonality condition w^v^ 7^ 0, known from ordi- 
nary matrix eigenvalue problems, to exponential matrix eigenvalue problems of the type 
K{?i, jLt) V = 0. The proof of Theorem 3.2 is entirely based on the standard calculus argu- 
ments for branching of solutions to algebraic systems as can be found in textbooks [1]. 
The details of the proof can be found in Section 6. 

The statements in Theorem 3.2 should be compared to two previous works considering 
the same situation: the branching of periodic orbits from an equilibrium losing its stability 
in FDEs with state-dependent delays. The PhD thesis of Eichmann (2006, [5]) proves the 
existence of the curve 13 >— > (x(^), pi(^), coil3)) and that it is once continuously differen- 
tiable (assuming only EC^ smoothness of/). Since jLt^(O) = due to rotational symmetry 
(see proof in Section 6) this is not enough to determine if the non-trivial periodic solutions 
exist for pi > or for jj, < for small (3 (the so-called criticalitry of the Hopf bifurcation, 
which is of interest in applications [15, 16]). Moreover, the non-resonance condition in 
[5] is slightly too strong, requiring that icog is the only purely imaginary root of detK{X, 0) 
(this assumption is different in the summary given in the review by Hartung et al. [12]; 
note that the publicly available version of [5] has a typo in the corresponding assumption 
LI), and only the pure-delay case (where the time interval J equals [ — t, 0]) was consid- 
ered. However, the techniques employed in [5], based on the Fredholm alternative, are 
likely to yield exactly the same result as stated in Theorem 3.2 if one assumes general EC'^ 
smoothness with k>2 (the formulation of EC'^ smoothness is already rather convoluted in 
[5]). 

Hu and Wu [13] use S^-degree theory [7, 18] to prove the existence of a branch of 
non-trivial periodic solutions near (x, jj,, co) — (0, 0, cOq). This type of topological methods 
gives generally weaker results concerning the local uniqueness of branches of periodic 
orbits or their regularity, but they require only weaker assumptions ([13] still needs to 
assume EC^ smoothness, though). Degree methods also give global existence results by 
placing restrictions on the number of branches that can occur. 
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4. Conclusion, applications and generalizations 



Periodic boundary-value problems for functional differential equations (FDEs) are equiva- 
lent to systems of smooth algebraic equations if the functional / defining the right-hand 
side of the boundary- value problem satisfies natural smoothness assumptions. These as- 
sumptions are identical to those imposed in the review by Hartung et al. [12] and do not 
exclude FDEs with state-dependent delay. There are several immediate extensions of the 
results presented in this paper. The list below indicates some of them. 

Further potential applications of the Equivalence Theorem 2.5 

Theorem 3.2 on the Hopf bifurcation is not the central result of the paper, even though 
it is a moderate extension of the theorem proved in [5]. Rather, it is a demonstration 
of the use of the Equivalence Theorem 2.5. The main strength of the equivalence result 
stated in Theorem 2.5 is that it permits the straightforward application of continuation 
and Lyapunov- Schmidt reduction techniques to FDE problems involving periodic orbits of 
finite period, regardless if the delay is state dependent, or if the equation is of so-called 
mixed type (that is, positive and negative delays are present). A source of complexity, for 
example, in [5, 7, 10, 13, 18, 26], is that techniques such as Lyapunov-Schmidt reduction 
or S^-degree theory had to be applied in Banach spaces. Theorem 2.5 removes the need 
for this, reducing the analysis of periodic orbits to root-finding in M"'^'-^'^^^-'. For example, 
Humphries et al. [14] study periodic orbits in FDEs with two state-dependent delays 
numerically using DDE-Biftool [6], alluding to theoretical results about bifurcations of 
periodic orbits that have been proven only for constant delay. Fist of all, Humphries et al. 
[14] continue branches of periodic orbits. Theorem 2.5 makes clear when these branches as 
curves of points (x, co, ji) in the extended space C°(T; E") x E x E are smooth to arbitrary 
degree: the Jacobian of (22) with respect to (py, co, /i) along the curve has to have full rank. 
Along these branches Humphries et al. [14] encounter degeneracies of the linearization 
and conjecture the existence of the corresponding bifurcations (backed up with numerical 
evidence) such as: fold bifurcations, period-doubling bifurcations, or the branching off of 
resonance surfaces (Arnol'd tongues) when resonant Floquet multipliers of the linearized 
equation cross the unit circle [20] . The Equivalence Theorem 2.5 provides a straightforward 
route to proofs that these scenarios occur as expected. Similarly, Theorem 2.5 will likely 
not only simplify the proofs about bifurcations of symmetric periodic orbits such as those of 
Wu [26], but also extend them to the case of FDEs with state-dependent delays. As long as 
one considers branching of periodic orbits with finite periods, the problem can be reduced 
locally (and often in every ball of finite size) to a finite-dimensional root-finding problem. 
This transfers also a list of results of sjonmetric bifurcation theory found in textbooks [8] 
to FDEs with state-dependent delays. 

Globally valid algebraic system 

The main result was formulated locally in the neighborhood of a given Xq e C^'^(T; E") 
and required only local Lipschitz continuity. The proof makes obvious that the domain of 
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definition for the map X, which maps between the function space and the finite-dimensional 
space is hmited by the size of the neighborhood of Xq for which one can find a uniform 
(£C) Lipschitz constant of the right-hand side F. In problems with delay the right-hand 
side is typically a combination of Nemytskii operators generated by smooth functions and 
the evaluation operator ev : C°(T; E") x T ■-^ E", given by ev(x, t) = These typically 
satisfy a (semi-) global Lipschitz condition (see also condition (Lb) in [12]): for all R there 
exists a constant K such that 

\fix)-fiy)\<K\\x-y\\o 

for all X and y satisfying \\x\\i < R and \\y\\i < R. Under this condition one can choose for 
any bounded ball an algebraic system that is equivalent to the periodic boundary-value 
problem in this bounded ball. For periodic orbits of autonomous systems this means that 
one can find an algebraic system such that all periodic orbits of amplitude less than R and 
of period and frequency at most R are given by the roots of the algebraic system. 



Implicitly given delays 

In practical applications the delay is sometimes given implicitly, for example, the position 
control problem considered in [23] and the cutting problem in [15, 16] contain a separate 
algebraic equation, which defines the delay implicitly. In simple cases these problems can 
be reduced to explicit differential equations using the standard Baumgarte reduction [2] 
for index-1 differential algebraic equations. For example, in the cutting problem the delay 
T depends on the current position x via the implicit linear equation 

Tit) ^ a - bxit) - bxit - Tit)), (28) 

which can be transformed into a differential equation by differentiation with respect to 
time: 

. ^ ^ _ -bvit) - [Tit) - a + bxit) + bxit - Tit))] 

1 + bvit -Tit)) ^^^^ 

(v(t) = x(t) is explicitly present as a variable in the cutting model, which is a second-order 
differential equation). The original model accompanied with the differential equation 
(29) instead of the algebraic equation (28) fits into the conditions of the Equivalence 
Theorem 2.5. The regularity statement of the Equivalence Theorem guarantees that the 
resulting periodic solutions have Lipschitz continuous derivatives with respect to time. 
This implies that the defect d = t — (a — bx — bx(t — t)) occurring in the algebraic 
condition (28) satisfies d(t) = — d(f ) along solutions. Since the solutions are periodic the 
defect d is periodic, too, and, hence, d is identically zero. The denominator appearing in 
Equation (29) becomes zero exactly in those points in which the implicit condition (28) 
cannot be solved for the delay t with a regular derivative. 

The same argument can be applied to the position control problem as long as the object, 
at position x, does not hit the base at position — w (the model contains the term |x -I- w| in 
the right-hand side) . 
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Neutral equations 



The index reduction works only if the delay t, which is itself a function of time, is not 
evaluated at different time points. For example, changing bx{t — T(t)) to bx{t — T{t — 1)) 
on the right-hand-side of (28) would make the index reduction impossible. However, 
certain simple neutral equations permit a similar reduction directly on the function space 
level. Consider 

- [A,(x + g(x))] =/(A,x) (30) 

where the functional / satisfies the local EC Lipschitz condition, defined in Definition 2.2, 
in a neighborhood [/ of a point Xq e C^'^(T;E"), and g : C°(T;E") ^ E" has a global 
(classical) Lipschitz constant less than unity (this excludes state-dependent delays in the 
essential part of the neutral equation). Then one can define the mapXg(j) as the unique 
solution X of the fixed point problem 

xit) = y(f) - g(AtX) for y near Jq = Xq + g(xo), 

which reduces (30) to the equation 

m = /(A,X,(y)) = /(X,(A,y)). (31) 

Equation (31) satisfies the conditions of the Equivalence Theorem 2.5. One implication of 
this reduction is that periodic solutions of (30) are k times continuously differentiable if 
the functional / is EC^ smooth in the sense of Definition 2.1 and g isk times continuously 
differentiable as a map from C°(T; M") into E". 



5. Proof of the Equivalence Theorem 2.5 

Theorem 2.5 is proved in three steps. First, we establish the existence of a locally unique 
fixed point of the map M^(-,p) using Banach's contraction mapping principle. This step 
requires only local EC Lipschitz continuity in the sense of Definition 2.2. In the second 
step we prove continuous differentiability of the map X and the right-hand side g of the 
algebraic system assuming that / is EC^ smooth. In the final step we prove higher-order 
differentiability, assuming that / is EC^ smooth for degrees k up to j^^. 

5.1. Decay of Fourier coefficients for integrals and smooth functions 

The following preparatory lemma states the well-known fact that, colloquially speaking, 
integrating a function makes its high-frequency Fourier coefficients smaller. In the fixed- 
point equation (18) of Theorem 2.5 the term Qj^L occurs, and we need this term to be small 
for large N. Recall that Qj^ removes the first N Fourier modes from a periodic function and 
Lx is the anti-derivative of x (after subtracting the average of x), see Equation (15) for 
the precise definitions. 
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Lemma 5.1 (Decay of Fourier coefficients of integrals) 

The norm of the linear operator Q^^L, mapping the space C^ (T; M") back into itself, 
is bounded by 

I log AT 
\\QmL\\,<C^ 

where C is a constant. The same holds in the Lipschitz norm (with the same 
constant C): 

logN 

llQiv-f-lki<c^. 



Proof We find the norm ||Qj^,L||o first, and start out with the well-known estimate for in- 
terpolating trigonometric polynomials for continuous functions on T. Let x be a continuous 
function on T with modulus of continuity co. Then (see [17]) 

llQ^xllo < CoCO j logN 

where Q is a constant that does not depend on x or N. A function co : [0, oo) ^ [0, oo) is 
called a modulus of continuity of a continuous function x if 

\xit)-xis)\<coi\t-s\). 

holds for all s and t e T. For a function x e C°(T; E") the anti-derivative 



[Lx](0 



x(5) —RqX ds 



has the Lipschitz constant ||x||o = max{|x(t)| : t e T} such that a modulus of continuity 
for Lx is co{h) — WxWqH. Consequently, 

||Q^Lx||o<Co^^^||^logiV, (32) 

where Q does not depend on x or N. This proves the claim of the lemma for j — 0. For 
X e C-'(T;]R") we notice that all derivatives of x up to order j are continuous. Applying 
estimate (32) to each of the derivatives of x we get 

||Q^Lx«||o < ^ logiV ||xW||o for l = 0,...j. 

Consequently, the maximum of the left-hand sides over all I G {0, . . . ,j} must be less than 
the maximum of the right-hand sides: 

2nCn logN 
WQ^LxWj = max ||Q^Lx«||o < logN max ||x«||o = 27iCo — ||x||^.. 
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which imphes the desired estimate for ||QjvL||j using the constant C = 2nCo. 

The estimate of Qj^L in the Lipschitz norm is a continuity argument. The operator Q^L is 
bounded (and, thus, continuous) on C°'-^(T; E"). For every element y of C^(T; R") (which is 
a dense subspace of C°'^(T; E")) the Lipschitz constant is identical to ||y'||o = max^gj |y'(f )|, 
and, thus, ||y||i = ||y||o,i- Let x„ e CHT;E") be a sequence of continuously differentiable 
functions that converges to x e C°'^(T; M") in the || • ||o ^-norm: ||x„ — x||o ^ ^ for n ^ oo. 
Then 

logN log AT 

llQw-f'-^nllo,! — WQn^^iiWi — ll-^nlll — ll-^nllo,l- 

On both sides of the inequality the limit for n — > oo exists, resulting in the desired estimate 
for llQjvLllo,!- □ 
A direct consequence of Lemma 5.1 is that the Lipschitz norm of Qj^x, ||Qw^llo,i5 go^s to 
zero for ^ oo for elements of C^'^(T; M"), so, for example, for a solution x of a periodic 
BVP: 

logN , logN 
llQiv^llo,! = llQivi^ llo.i < C^llx'llo,! < C^llxlh 1. (33) 



N " N 



5.2. Proof of Splitting Lemma 2.4 



For any given integer N > we have to show that the pair (x,p) e C°(T; E") x ]R"x(2iv+i) 
satisfies 

r 

x(t) = x(0) + F(x)(5) d5 for all t e T, and (34) 
Jo 

p = R[^x, (or, equivalently Ej^p = Pf^x) (35) 

if and only if it satisfies the system 

X ^ E,^p +Qf^LF{x), (36) 
= R^ [PoFix) + Qo[E^p-Pj,LFix))], (37) 

Assume that x e C°(T;E") satisfies (34), and let p =Ri^x. Subtracting equation (34) 
for t = — 71 from (34) for t — n implies that the average of F(x) is zero. Thus, RqF{x) — 
and PqF(x) — 0. Since Ly — y(s) —R^yds, the identity (34) implies (in combination 
withRof (^) = 0) 

x(t) = x(0) + [LF(x)](f). (38) 

Applying projection Qjv to this identity we obtain QjyX = QjvLF(x). Adding (35) to 
this we obtain equation (36). Applying projection PjvQo (which is the same as Qo^'w) to 
(38) we obtain QqPn^ — Qo^w^^C-^)- Inserting E^^p for P^^x into this identity leads to 
Qo[FjvP — ■Pjv-f'f^C-^)] = 0- Since PqF{x) — 0, this in turn implies (37). 

Applying Pj^ to (36) implies P^x = E^^p (and p =Rf^x) immediately. The expres- 
sion inside the parentheses of R^^ in equation (37) is a sum of two parts that each have to 
be zero (since they are both in the image of Pjy on which Rf^ in injective). The projection 
Qo subtracts the average from its argument. Hence, (x,p) e C°(T;E") x ]R"x(2iv+i) satisfies 
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(36)-(37) if and only if there exists a constant c e E" such that the triple (c, x,p) satisfies 
the system of equations consisting of (36) and 

= RoFix) (39) 
EqC = E^p - Pf^LF(x). (40) 

Note that Eq maps the constant c e M" to a function that equals this constant for all t e T. 
In this system, (39) ensures that the average of F(x) is zero. Equation (40) is an equation 
in the finite-dimensional space rgP^^- Subtracting (40) from (36) gives 

X — EqC + LF{x). 

This equals (34), keeping in mind that [Ly](t) = j^jis) - RoJ ds ([Ly](0) = for all 
y e C°(T; M"), hence, x(0) = c), and using RoF^x) = (see equation (39)). □ 



5.3. Unique solvability of the fixed point problem (18) 

Let Xq be an element of C^'^(T; E"), for example, a solution of the periodic boundary value 
problem x(t) = /(A^x) = F(x)(t). Consider a closed ball B°'^(xo) of radius 5 around Xq 
in the Lipschitz norm: 

Bfix,) = {x e C°'HT;E") : ||x - Xollo,i < 5}. 

The superscript "0, 1" indicates which norm is used to measure the distance from Xq and 
that only elements of C°'HTr;M") are included. 

Lemma 2.3 implies that F is Lipschitz continuous with respect to the || • Hg-norm in 
B°'^(xo) if we choose 5 sufficiently small (thus, F is also called EC Lipschitz continuous in 

Bl-\xo)y. 

linx)-F(y)||o<K||x-y||o (41) 

for all X and y in B°'^(xo) and a fixed i^C > 0. In any ball in which F is EC Lipschitz 
continuous, F is also bounded in the Lipschitz norm: 

||F(x)||o,i</? forallxeB°'^(xo). (42) 

See Lemma A.5 in Appendix A for the proof. 
We can now formulate a lemma about the unique solvability of the fixed point problem 

X ^Ef^p+Q^LFix). 

This unique solvability and the Splitting Lemma 2.4 allow us to reduce the periodic BVP 

x(t) = /(AfX) to a system of algebraic equations. Remember that E^^p takes a vector p 

of 2N + 1 Fourier coefficients and maps it to the periodic function having these Fourier 

coefficients, R^x extracts the first 2N + 1 Fourier coefficients from a periodic function x, 

PffX projects the periodic function x onto the space spanned by the basis b_^, . . . , bjv 

Qf^ = / — Pjy sets the first Fourier modes of a function to zero. (P^ and Qj^ are projections 

in the function space, and Rf^ and ni^P between the finite-dimensional subspace rgPjv 
and E"^f2jv+i)_) 
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Lemma 5.2 (Unique solvability of fixed point problem) 

Let Xq be in C^'^(T; E"), and let 5 > be such that 

I ||F(x)||o,i<K and ||F(x) - F(y)||o < K||x - y ||o (43) 

for all X and y e 5°^ (xq) and for some constants K > and R > depending on 

15. Then for any sufficiently large AT the fixed point problem 
x = E^p + (l^h¥{x) (44) 

has a unique solution x e 5°^ (xg) for all vectors p e ]r"x(2n+i) ^j^g neighborhood 
U of KjvXq given by 

[/ = {p e E-x^^'^+i) : ||£^ [p -i^^xo] 11,^, < 25} . (45) 

Moreover, if x e B°'^(xo) is continuously differentiable and satisfies x' = F(x) then 
its projection p = RjyX is in the neighborhood U , and x and p satisfy (44). 

Note that U is an open set of M"'^'^^'^"'"^^ since is an isomorphism between rgP^^^ equipped 
with the II • llo i-norm, and ]R"'^'^2jv+i) ^Ve have to prove the unique solvability of the fixed- 
point problem in a slightly larger ball (radius 65) and for a slightly larger range of 
parameters p (note the 25 in (45)) in order to establish one-to-one correspondence in the 
ball of radius 5. 



Proof The idea is, of course, that the function 

M^(-,p):x^£^p+QjvLF(x) 

maps the closed ball B°^^(xo) back into itself and is uniformly contracting for suitably large 
N and vectors p gU. 

First, any closed ball B°'^(xo) is closed (and, thus, forms a complete metric space) 
with respect to the || • Hg-norm. This completeness is a simple continuity argument: let 
y„ = Xg +z„ be a Cauchy sequence in B°'^(xo) with respect to the || • ||o-norm. Then z„ 
converges to a continuous function z, and, since Wz^Wo < l|Zn||o,i < r, for all n, the maximum 
norm of z is also bounded by r: ||z||o < r. We only have to show that the Lipschitz constant 
of z is bounded by r, too. Let e > be arbitrary and let t 7^ 5 be arbitrary in T. We select 
some n such that ||z — z„||o < e|t — s|/2. Then 

|z(t) - Z(5)| < |z(f ) - Z,(t)| + Kit) - Z„(5)| + |Z,(5) - Z(5)| 

< £\t — s\ + r\t — s\ < (r + £)\t — s\. 

Thus, the Lipschitz constant of z is less than r + s for arbitrary e > 0. Hence, ||z||o 1 < r, 
completing the argument for completeness of B°'^(xo) with respect to the || • ||o-norm. This 
completeness implies that we can apply Banach's contraction mapping principle in a ball 
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B°'Xxo), a ball of Lipschitz continuous functions, using the weaker maximum norm in the 
following. 

We choose the radius r of the ball equal to 65 (5 was chosen in the lemma such that the 
estimates (43) are true for the constants K and R), Thus, 5°^ (xg) is the set to which we 
want to apply Banach's contraction mapping principle. To ensure that the map Mjv(-,p) 
maps into 5°^ (xq) for p & U, and that Mjv(-,p) is a contraction we pick N large enough. 
Specifically, we pick N such that 

25 

llQ^Xollo,! < 25, ||Q^L||o,i<— , 

1 logN , , , , ^46) 

llQ^Lllo < — , < l/max{l, + ||xolli,i) /5}, 

where R and K are the bounds on F given in the conditions of the lemma, in Equation (43). 
We know that these bounds exist due to Lemma 2.3 (see Equation (41)) and Lemma A.5 (see 
Equation (42)). We know that choosing N according to (46) is possible from Lemma 5.1 
and estimate (33) following Lemma 5.1. 
Let us check first that x Ef^p + Q[^LF{x) maps the closed ball B°g (xq) back into itself: 

\\EnP + QnLF(x) - Xo\\^^^ < 

< ||£jv [p-RnXo] -QjvXo + Qw-f.F(x)||g_^ 

< \\e^ [p -RnXo] + ||Q]v^o||o,i + ||Qiv-^||o.i ll-^WIIo.i 

25 

<25 + 25 + — R = 65. 
R 

Here we used the bounds (46) implied by our choice of N and the definition (45) of the set 
U of permitted p, and the bound on ||F(x)||o j, which is determined in (43) by our choice 
of 5. 

Second, let us check that x Ef^p + QjvLF(x) is a uniform contraction in 5°^ with 
respect to the || • ||o-norm: 

\\q^L [Fix) - F(j)] ||„ < IIq^L \l ||f(x) - F(y)\l < ^K\\x - j||o < ^||x - y ||o. 

Again, we exploited the bounds (46), implied by our choice of N, and the Lipschitz constant 
K oi F determined in (43) by our choice of 5. 

Since 5°^ (xq) is complete with respect to the || • ||-norm Banach's contraction mapping 
principle implies that the fixed point problem (44) has a unique solution x e B°'^(xo) for 
pGU. 

Finally let us check that for x e B°'^(xo)nCHT; M") satisfying the periodic BVP x' = F(x) 
the projection p — R^^x is in U. For this we have to prove that if ||x — XqHoi < 5 and 
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x' — F(x) then ||Pjv(x — Xo)||o,i < 25. We can estimate ||Pjv(x — Xo)||o,i via 

||P^(x - Xo)||o,i < ||(/ - Qiv)U - ^o)llo,i < Ik - ^ollo.i + llQivU - ^o)llo,i (47) 

logN 

<5 + C^||x-Xolli,i (48) 

< 5 + C-|^max{|x -Xollo,i,lk -^ollo.i^ (49) 

< 5 + C — max{5, ||x'||o,i + \\x%^,} (50) 

logN 

= 5 + C^max{5, ||F(x)||o,i + ||xolli,i} (51) 
logN 

<5 + C-^max{5,R + ||xo||i,i}<25. (52) 

The inequahty (47) follows from the definition of Pj^ and and the triangular inequality 
for the II • llo i-norm. The step from (47) to (48) uses the estimate (33) for the norm 
llQjvJIlo,! for elements y of C^'^(T;E"). It also bounds ||x — Xgllg ^ by the radius 5 of the 
ball. Step (49) splits up the || • ||^ j norm into its two parts which are estimated separately 
in the following steps. One part, ||x — XqIIoj is bounded by 5 (the radius of the ball), the 
difference of the derivatives is bounded by a triangular inequality for its parts, ||x^||o i and 
llxgllo 1 in (50). To get to (51) we use that x satisfies the BVP x' = F(x). We also bound the 
norm of x^ by llxglli i. Finally, in (52) we estimate the Lipschitz norm of F(x), ||F(x)||o^ by 
the bound R determined in (43) by our choice of 5. The right-hand side of (52) is (strictly) 
less than 25 by our choice of N, see (46). □ 



5.4. Lipschitz continuity of the algebraic system 

The Splitting Lemma 2.4 guarantees in combination with the unique existence of the 
fixed point of Mjv(-,p), proven in Lemma 5.2, the equivalence between the periodic BVP 
x(t) = /(AfX) and the algebraic equation g(p) = for x inside the ball ^(xg), where g 
is given in (19) by 

g:pGU^R^ [PoF(X(p)) + Qo (F^p - P^LF(X(p)))] G E"x(2iv+i)^ ^^ere (53) 
X :p e [/^ C°(T;E"), and X(p) is the fixed point of Mjv(-,p) in B°5^(xo). (54) 

The relation between p G U and x e B°'^(xo) is given via p = Rp^x and x = X(p): if x 
satisfies the periodic BVP then p —R^^x satisfies g(p) = 0, and, vice versa, ii p G U satisfies 
g(p) = then X = X{p) satisfies the periodic BVP The domain of definition, D(X) = U 
is an open set, however the map X (and, thus, g) can be extended continuously to the 
boundary of U: Mj^{-,p) maps into the closed ball 5°^ back into itself also for p on the 
boundary of U and it still has contraction rate 1/2 with respect to the || • Hp-norm. 

The remainder of the section addresses the remaining open claim of the Equivalence 
Theorem 2.5, namely the regularity of the maps X and g. Using only local EC Lipschitz 
continuity (Definition 2.2) we can prove the Lipschitz continuity of g and X: 
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Lemma 5.3 (Regularity of X and algebraic system) 

1. For all p in the neighborhood U = D{X), defined in (45), the image X{p) is 
in C^'^(T;E") (that is, X(p) e CHT;E") and its time derivative is Lipschitz 
continuous), 

2. X is Lipschitz continuous with respect to the || • ||i-norm for its images: there 
exists a constant Cjv such that 

\\Xip) -Xiq)\\, < C^\p-q\ for all p and q in U, 

3. the map p gU ^ [RoFiXip)), Pf^LFiXip))] e E" x ^"^(^iv+i) ypschitz 
continuous in U. 



Proof For a function y e rgP^^j differentiation is a bounded operator: y' = Dj^y. The 
vector Rj^y of the first 2N + 1 Fourier coefficients of a function y and the vector i?^ [y'] 
satisfy Rjv[y^] = Dj^Rj^y where D^^ is a matrix (independent of y). Hence, y' — Ejv^w-^jvJ 
for all y e rgPjv such that we can define Df^ — EjvDjv^jv- Denote X(p) as x. By definition 
of the map X, x — E^p + Qj^LF{x). The right-hand side of this fixed-point equation is 
differentiable with respect to time, giving 

x' ^Df,E^p+QoF(x)-Df,Pf,LF(x). (55) 

This guarantees that x e C^(T;]R"). Equation (42) ensures that ||F(x)||o,i < R, which 
implies that the right-hand side of (55) is Lipschitz continuous in time. This in turn implies 
that x' is Lipschitz continuous in time (thus, x e C^'^(T;E")), and 

||x1|o,i < |K£jvllo.ilPl + llQollo,i^ + l|i^ivPivi'llo.i-R- 

Representation (55) also implies point 2: let x =X{p) and y —X{q) be two functions in 
the image of X: 

||x' - y'llo < ||D^£ivllolP -q\ + (llQollo + \\D^PnL\\o)K\\x - y||o, (56) 

where K was the EC Lipschitz constant of F in B^'^{xq). The difference x — y in the 
II • llo-norm is bounded due to the contractivity of the right-hand side in fixed point problem 
(44) defining X (the || • ||o-norm was the metric used to apply the contraction mapping 
principle) : 

1 

||x -yllo < llFwIlolp -ql + llQw-f'[F(x) -F(y)||o < ||£;vllolP - qI + 2!!^ ~3^llo- 

Thus, 

||x-y||o<2||£jvllolp-qL 
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which, combined with (56), gives Lipschitz continuity of X as a map from U into C^(T; M"): 
||x - y'llo < [\\DM\o + (llQollo + \\D^PnL\\o)2K\\EM \p - q\ =: Q|p - q\. (57) 

Point 3 is a direct consequence of the Lipschitz continuity of F with respect to || • Hg-norm 
in (xq), the Lipschitz continuity of X on U in the || • Hg-norm, and the fact thatX maps 
intoB°5^(xo). □ 



5.5. First-order differentiability of the algebraic system 

Until now we have only used the EC Lipschitz continuity (in the sense of Definition 2.2) of 
the right-hand side F in the ball (xg) with respect to the || • ||o-norm. We can expect 
that the right-hand side g of the algebraic system, defined in (53), is smooth only if we 
require more smoothness of the right-hand side / (which enters F in the algebraic system). 

We first discuss first-order differentiability of the map X and the right-hand side g, de- 
fined in (53) and (54). For this we assume EC^ smoothness of f as defined in Definition 2.1. 
For X e C^T; E") n Bgg (xq) the norm of fix, •) as an element of L(C°(T; E"); E") (the 
space of linear functionals mapping C°(T; E") into E") is less than or equal to K, the EC 
Lipschitz constant of F (and, hence, /) in (xg) assumed to exist in the conditions of 
Lemma 5.2. 

Let us define the map 

d^F :C\T;E") x C°(T;E") ^ C°(T;E"), [a^F(v,w)] (t) = d^fiA,v,A,w). 

If V e C^(T;E") and w e C°'HTr;E") then the map d^F defined above is indeed the 
derivative of F in v with respect to the deviation w (see Lemma A. 6 in Appendix A): 

||F(v + w)-F(v)-aiF(v,w)||o ^ 
lim — = 0. (58) 

Ikllo.i^o 

Part of the definition of EC^ smoothness for / is that the map d^f is continuous in both 
arguments, v e C\T; E") and w e C°(T; E"). One can then apply Lemma A.4 to d^f to 
conclude that the map d^F (a composition of and d^f) is continuous with respect to the 
II • llo-norm in its image space as a map of both arguments (in their respective norm), too. For 
V e (xq) the norm of the linear map 5^F(v, •) as an element of L(C°(T; E"); C°(T; E")), 
the space of continuous linear functionals from C°(T; E") back to itself, is bounded by the 
EC Lipschitz constant K of F in Bg^ (xg). 

The additional regularity assumption on / and its implications for F permit us to improve 
our statements about regularity ofX and the algebraic system: 

Lemma 5.4 (Continuous differentiability of X and the algebraic system) 

Assume that the right-hand side / is EC^ smooth in the sense of Definition 2.1. 
Then the regularity statements about the map X, defined in (54), and the right- 
hand side of the algebraic system, defined in (53), can be extended: 
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1. X(p) is in C^(T;]R") for a\\pGU = D(X), the domain of definition of X, and 
p X{p) is continuous with respect to the || • ||2-norm for its images. 

2. The map X, which maps U into C^(T; E") according to Lemma 5.3, is contin- 
uously differentiable with respect to its argument p using the || • ||^-norm for 
its images. 

3. The map p gU ^ [RQF{X{p)),Pt^LF(Xip))] e M" x ^"^^2^^+^^ is continu- 
ously differentiable with respect to p. 



Proof Let p e [/ = D(X) c ^"^^2^+^^, where U is defined in (45), and let us denote X(p) 
by X. Lemma 5.3 ensures already that x is in C^'^(T;]R"). Lemma A.7 in Appendix A 
proves that F(x) e C\T; E") for x e C^{T; E") (choosing D = C°(T; M") and = in the 
assumptions of Lemma A.7). This implies the first statement, thatX(p) e C^(T;E"): since 



andX(p) e C^'^TjE") (see Lemma 5.3), F(X(p)) is in C^(T;M"), and, thus, LF(X(p)) 
is in C^(T;E"). Hence, X(p) is an element of C^(T;E"), too. Furthermore, Lemma A.7 
states that F is continuous as a map from C^(T; E") into C^(T; M"). Since X is continuous 
as a map from U into C^(T;E") (in fact, it is Lipschitz continuous, see Lemma 5.3), the 
right-hand side of (59) inp is continuous with respect to the || • H^-norm. This proves the 
first point. 

Concerning the second statement: again, let po be in 17 = D{X), and choose a small open 
neighborhood f/(Po) which has a positive distance to the boundary of U. We will prove 
point two for all p e f/(Po)- Let us choose an initial £q sufficiently small such that p + hq 
is still in U for h e {—Sq, Sq), all p e U{pq), and all q with \q\ < 1. Let us introduce the 
difference quotient for h e {—Sq, Sq) \ {0}: 



The maps z maps [(-^o, ^o) \ {0}] x 5^(0) x [/(po) c E x E"x(2iv+i) ^ ^nx(2jv+i) ^^^^ 
C\T; E"). We first prove that z has a limit for h ^ in C^(J; E"), and that this limit is 
achieved uniformly for all p e f/(Po) and |q| < 1. By definition of X, z satisfies the fixed 
point equation (dropping all arguments from z) 



Xip)^E^p+Q^LFiXip)) 



(59) 



z(h,q,p) = -[Xip+hq)-Xip)]. 



(60) 




(61) 



for h e (— ^o) \ {0}. Let us introduce 



Ai(p,z,h) = 



i[F(X(p) + hz)-F(X(p))] ifh/0 
aiF(X(p),z) ifh = 0. 



(62) 
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which maps [/(po) x C°'HT;E") x R into C°(T;E"). The hmit (58) imphes that A i is 
continuous in all arguments (insert v — x,w — hz into (58)). UsingAj we extend the fixed 
point problem (61) to h = 0: 

z = £^q + QjvIAi(X(p),z,h). (63) 

The following intermediate lemma proves that the fixed point problem (63) has a unique 
solution: 

Lemma 5.5 (Fixed point problem for linearization) 

There exists an e > and constants Cq> and > such that the map 

y.zG C°'HT; E") ^ + Q^LA^iXip),z, h) e C^'^T; E"), 

which depends on the additional parameters p, q and h, has a unique fixed point 
z* in 

B = {z e C°'HT;M") : ||z||o < Cq and ||z||o,i < Cj 
for all h e (-£, s), all p e Uipo) c [/ = D(X) c m"x^2]v+i) ^jj ^ ^ ^n(2]v+i) ^j^j^ 
Iql < 1. The fixed point z^ is an element of C^(T;]R") and depends continuously 
on h, p and q with respect to the || • H^-norm. 

Note that the s we have to choose in Lemma 5.5 is smaller than the initial £q for which the 
difference quotient z is defined. 

Proof of Lemma 5.5 First of all, since is continuous in all arguments, the map y is 
continuous. Moreover, since x' = (X{p)y and x = X{p) depend continuously on p (see 
Lemma 5.3 and expression (55)), the map y also depends continuously on the parameters p, 
q and h (that is, the expression E„q + Q^LAi{X{p),z,h), defining y, depends continuously 
on z, p, q and h with respect to the || • ||o ^-norm). We choose the constants Cq > and 
Ci> such that 

Co>2||£jvllo (64) 
+ (llQollo + l|i?iv^ivi'llo)^^Co, (65) 

where K is the Lipschitz constant of F with respect to the || • Hg-norm in 5°^ . 

We choose £ < Sq such that for all z satisfying ||z||o i < and all p e f/(Po) the function 
X(p) + hz is in 5°^ (xq) for all h e (— e, e). This implies that for any z^ and Z2 satisfying 
1 1^1 1 1 0,1 — Q and ||z2llo,i < Q we have 

UwiXip) + hz,) - F(X(p) + feJIlo < i^llzi - Z2II0, 

h (66) 

\\d'FiXip),z,-Z2)\\o<K\\z,-Z2\\o, 
where K was the Lipschitz constant for F in (xg), and, thus, 

1 

llr(zi) - r(z2)llo ^ 2"^! ~^2llo (67) 
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for all h e (—s, s) by choice of N {N was such that ||Q^L||o < (,2K) ^). This estimate for y 
implies 

llr(z)llo < ll^ivllo + ^llzllo if l|z|lo,i < Q, (68) 

since ^(O) = E^^q and |q| < 1. Moreover, the two inequalities (66) imply that for h e {—£, e), 
||z|lo,i ^ Ci and p G U{pq) the maximum norm oiAi{p,z,h) is bounded by iiC||z||o: 

\\A,(p,z,h)\\o<K\\z\\o (69) 

The time derivative of ^(z) exists and its || • ||o-norm can be estimated by differentiating the 
expression E^q + Q}^LA^(X(p),z,h), defining y, with respect to time in the same manner 
as we obtained (56) (we insert (69) to bound ||Ai(p,z,h)||o): 



d 



< llD^E^IIo + (llQollo + \\Df,P^L\\,)K\\z\\o. (70) 





The combination of the bounds (68) and (70) and the definition of the constants Cq and 
Ci guarantee that 'y(z) maps the set 

B = {z e C°'HT;M") : ||z||o < Cq and ||z||o,i < 

back into itself. The contraction estimate (67) for the || • ||o-norm and the completeness of 
B with respect to the || • Hg-norm make the contraction mapping principle applicable with a 
uniform contraction rate for all p e U{pq), all |q| < 1 and h e {—£, s). This ensures that 
the fixed point z^ depends continuously on p, q e ]^"(2w+i) ^ ^ {—£, e) with respect to 
the II • llo-norm (since y is continuous with respect to z, h, q and p). 

The time derivative z^ of z^ also exists and is continuous in p, q and h: we differentiate 
the fixed point equation (63) for z^ with respect to time (in the same way as done in (55)) 
to get 

zl =D^E^q+QoA,iXip),z„h) - D^P^LA^iXip),z„h), (71) 

which is a continuous function in p, q and h with respect to the || • ||o-norm (note that 
z^ depends on p e U{pq), q and h). Thus, the fixed point z^ is in C^(T; E") and depends 
continuously on p, q and h with respect to the || • ||i-norm. □ 

Proof of Lemma 5.4 continued As a consequence of Lemma 5.5 we may write the fixed 
point z^ of 7 as a function of h, q and p: z^(h, q,p) maps h e {—£, s), q in the unit ball 
]^nx(2jv+i) p g Uipo) continuously into CHT;E"). It is also identical to z(h,q,p), 
defined in (60) as the directional difference quotient ofX. Thus, the directional difference 
quotient z(h,q,p) has a limit for h ^ in the || • H^-norm, and this limit equals z^(0, q,p). 
Moreover, this limit z^(0,q,p) depends continuously on p and q in the || • ||^-norm (as 
proved in Lemma 5.5), and it is linear in q (since A^(p,z, 0) is linear in z). Thus, z^(0, q,p) 
is the Frechet derivative: 

||X(p + q)-X(p)-z,(0,q,p)||i 
lim — = 0. (72) 

q-O |q| 



31 



Consequently, the map (p,q) >— > z*(0, q,p) = d^X{p)q is continuous in the || • ||i-norm as 
claimed in the lemma. 

The third statement of Lemma 5.4 is a consequence of the second statement and the 
fact that the difference quotient of F has a limit in the || • ||o-norm if it is taken between 
arguments in CHT;M") (see (58)). We split the difference quotients into two parts: 

FiXip + hq))-F{X{p)) FiXip) + hd'Xip)q)-FiXip)) , 



^ FiXip+hq))-FiXip) + hd'Xip)q) 



(73) 
(74) 



The right-hand side in (73) converges in the || • ||o-norm to d^F{X{p),d^X{p)q) for h — > 0, 
since X(p) and d^X{p)q are in C^(T;E") because F is EC^ continuous (see the second 
point of the lemma for the regularity of d^X(p)q and Lemma 5.3 for the regularity of 
X{p)). For the term in (74) we can apply the local EC Lipschitz continuity (all arguments 



are in B^'^Xq) for p e U{pq), |q| < 1 and h e (e, e)) such that we get 



F(X(p + hq)) - FiXip) + hdX{p)q) 



<K 



Xip + hq)-Xip) 



dXip)q 



which converges to for h — > due to the second statement of the lemma (K is the EC 
Lipschitz constant of F in 5°^ (xq)). Consequently, we obtain from the limit of (73) for h — » 
that the directional derivative of F(X(p)) in direction q is equal to d^F{X{p), d^X{p)q), 
which is continuous with respect to p and q and linear in q. Thus, 



d 



F(X(p)) 



q^d'FiXip),dXip)q), 



(75) 



and p F(X(p)) is continuously differentiable with respect to p in the || • ||o-norm. Note 
that we use the notation not enclosing q in the bracket in (75) to highlight that this is a 
classical derivative with respect to a finite-dimensional variable. The linear operators Rq 
and Pf^L preserve the continuity (and the linearity in q) of (75). □ 

Lemma 5.6 (Spectral radius of the linearized fixed point problem) 

For X =X(p) (where p gU — D(X)) consider the linear map 

M :z^Qjvl5^F(x,z). 

The spectral radius of M as a map from C°(T; M") back into itself, or as a map from 
CHT;M") back into itself, is less or equal 1/2. 



Proof Since M is compact as an element of L(C\T; W); C\T; E")) (the space of linear 
functionals from C'^(T;M") back to itself) for = and k = 1, the spectral radius is 
identical to the modulus of the maximal (in modulus) eigenvalue, which is of finite 
algebraic multiplicity if it is different from zero. An eigenvector z corresponding to this 
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maximal eigenvalue is an element of C^(T; E") such that the spectral radius of M is the 
same for k — and k—1. 

Since x and z are both in C^(T; E") we have that 

1 

5 ^F(x)z =lim- [F(x + hz)-F(x)] (76) 

For X — X{p) where p &U, and h sufficiently small the arguments of F, x + hz and x, both 
lie inside 5°^ such that the EC Lipschitz constant K applies to the difference: 

d|F(x + fe)-F(x)||o<K||z||o. 
h 

Since ||Qwi|lo < l/(2if), (76) and (5.5) combine to 

1 

\\Mz\\o < -llzllo- 

As z is an eigenvector corresponding to the largest eigenvalue, the spectral radius of M is 
less or equal 1/2. □ 
Thus, the derivative z — dX{p)q ofX in p is the unique solution of the contractive linear 
fixed point problem in C^iT; E") 

z^E^q+Q^Ld^FiXip),z). (77) 



5.6. Higher degrees of smoothness 

We observe that (x,y) = {X{p),d^X(p)q) satisfies the system of equations 

X ^Ef^p + Qf^LF(x) ^^^^ 

y ^Et^q+Qj^Ld'^Fix,y). 

This has a similar structure to the original fixed point problem (44) but in dimension 
= 2n with the variables (x,y) and parameters (p,q). Thus, we aim to apply a linear 
version of the arguments of Section 5.5 recursively, assuming that / is FC*^ smooth as 
recursively defined in Definition 2.1. Throughout this section we assume that / is EC^ 
smooth for all degrees up to order j^^^. 

For higher-order derivatives, we introduce the spaces Dj and the operators d^F for j > 
recursively: 

Do = C°(T;M") Dj^D]_^xDj_, 

d^F -.Dj ^ C°(T;M"), [a^F(x)](t) = d^fiA.x). 

The spaces Dj are products of the type (9), and the argument x of d^F and d^f is in Dj, a 
product of 2^ spaces. We also recall that the notion of subspaces D| of higher-oder (fc > 0) 
differentiability for product spaces such as Dj was introduced in Section 2. For example, 

= C''(T;E"), 

dI = D^+^ X = C^+\T; R") x C''+\T; E") x C^+\T; E") x C^T; R"), etc.. 
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all with their natural maximum norms. The maps d^F are all continuous and map indeed 
into C°(T; E") due to the continuity of d^f and (applying Lemma A.4 to Dj, d^F and 
d^f). It is also clear from the definition that d^'^'^F - d^[d''F] if j + < j^^^. We will also 
use the notation L(D|; Df) for the space of linear bounded functionals mapping from 

into D^. 

The following lemma is a consequence of the EC^ smoothness of/. 
Lemma 5.7 

^ For j + k< jinax the operator d^F is a continuous map from into C'^(T; E"). 

Proof of Lemma We have to apply Lemma A. 7 from Appendix A inductively over the order 
of differentiability (fc). To start the induction for = we can apply Lemma A.4 to Dj, d^F 
and d^f. For the inductive step let us assume that for k we know that d^F : D*^ C^(T; E") 
is continuous for all j < j^^^ — k. Let us fix a j < j^^^ — k — 1. We have to show that d^F 
maps D^"*"^ continuously into C^~^^(T; E"). We know (by inductive assumption) that d^F 
maps D| continuously into C*'(T;E") and that d^+'^F maps dJ^^ = dJ+^ x dJ continuously 
into C'^(T; E"). Thus, we can apply Lemma A.7 to d^F (this takes the place of the operator 
F in Lemma A.7) and D = D'^, obtaining that d^F : D|+^ C'^+HT; E") is continuous. □ 
An immediate consequence of Lemma 5.7 is thatX(p) and dX{p)q, as constructed in 
Section 5.5, are as smooth as the right-hand-side: 

Lemma 5.8 (Smoothness of X and dX in time) 

Let / be EC^'"'"' smooth. For every p gU = D{X) and every q e /?"(2jv+i) ^.j^^ func- 
tions X(p) and dX(ip)q satisfy X(p) e C^--+i(T;E") and dXip)q e C^--(T;E"). 
Moreover, the maps 

p^Xip) e C^--+HT;E") and [p,q] ^ dXip)q G C^">-(T;E") 

are continuous. 

Proof The function x = X(p) satisfies x = E^p + QjvIF(x). Since F maps = C^(T; E") 
back into itself continuously for all k < j^^^^, Qj^L maps into D^'^^ continuously for 
all k, and E^^p e C°°(T;E"), the fixed point equation implies the following: if x e Dg 
then F(x) e D^, thus, x = £jvP + Qw^-^^C-^) ^ I>o^nfor all < j^ax). Similarly z = 
Fjv<3 + Qw^^^f (x)z, and maps into for all k < j^^^ — 1. Thus, the fixed point 
equation implies: if z e and x e Dq"*"^ then (x,z) e D^, thus, 5^F(x,z) e D^, thus, 
z = EjvQ + QjvL5^F(x,z) e Dq"*"^ for all fc < — 1. All of the above dependencies are 
continuous such that the continuous dependence on p and q in the norms of Dq""^^^ and 
Dq""™, respectively, follows. □ 
We plan to find the derivatives of the map X inductively through fixed point equations 
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of the form (78) . In order to set up the recursion we define inductively the operators Fj by 
Fo(x) = F(x) for X e Do (79) 

^ eD, =Dj_, xD,. (80) 









d'Fj_,ix,y)_ 



for 



eD, =D' X D;. 

; j-i } 



Note that Fj is always linear in its second argument, y, since d^Fj_i is linear in its second 
argument. The operators Fj are combinations of derivatives of F. The plan is to study 
fixed-point problems of the type x = E^^p + Qf^LFj{x) (with j = 1 we obtain (78)). Before 
doing so, we establish which spaces the operators Fj map into: 

Lemma 5.9 (Image of right-hand side) 

For j + l + k < jujax the operator maps , continuously into D^. In particular, 
Fj maps Dj continuously back into itself. 



Proof The statement of the lemma follows inductively from the definition of Fj and D^. 
We apply Lemma 5.7 to start our induction over j (for j = the statement is identical to 
Lemma 5.7). For the inductive step let us assume that we know that d^Fj_i maps D^_,_;_i 
continuously into D'^_^ for all k and I satisfying l + k< — j + 1. By definition (80) of 
Fj the derivative d^Fj for I < j^^^ — j is 



d^F, 



y 



d'Fj_,{x) 
d'^'Fj_,{x,y) 



for 



D 



i+j 



D 



i+j-i 



X D 



i+j-i- 



The first component, d^Fj_i maps Df^^_^ continuously into D'^^I for all k from to jniax~'~i 
(this is the assumption of the inductive step when one shifts the index k down by 1). 



Similarly, 5'+^Fj_i maps Dj+;_^ x = D'^_^^ continuously into Dj_^ for all k from to 

Jmax ~l ~ j> again due to the assumption of the inductive step. Consequently, 5'Fj maps 



D 



k+1 
j+l-1 



X D 



j^i_^ = Dj_^^ continuously into Dj for all k from to j,^^^ — I — j, which is the 
statement we had to prove for the inductive step. □ 
Even though the map x e D j d^Fj^x, •) e L{Pj;Dj) is in general not continuous, the 
following map is: 

Lemma 5.10 (Continuity in operator norm) 

I For i < jmax the map x e Dj — > Qj^Ld^Fj{x,-) e L(D|;D|) is continuous with 
respect to x e D\ 



Proof of Lemma 5.10 The EC^ smoothness of / (for k < j^^^) implies that Fj is contin- 
uously differentiable (in the classical sense) as a map from Dj into Dj. Thus, the map 
X >-> d^Fj{x, •) as a map from Dj into L(Dj; Dj) is continuous. Recall that the operator L 
involves taking the anti-derivative of its argument such that L : Dj ^ D^. Since Q^^L maps 
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Dj continuously into Dj, the map x Q^Ld^Fj(x, •) is continuous as a map from D| into 
L(D];D]). ' □ 

The following theorem provides continuous differentiability of order j^^^ for X and the 
map p ^ F{X(p)) if the right-hand side is EC'^ smooth in the sense of Definition 2.1 for 

^ — imax- 

Theorem 5.11 (Smoothness of algebraic system and X) 

Define Uq — n{2N + 1) and rij — 2^ rig, and the maps 

Xo:pGU = DiX) c E"o ^x(p) e Do and 
7o : p e f/ = D(X) c E"o ^ F(X(p)) e Dq, 

and assume that / : Dg = C°(T;E") E" is EC-'""- smooth. Then the following 
maps exist and are continuous for all j up to j^^^: 

Xj : [p,q] G DiXj) := DiXj_,) x E".- C E". ^ [X^._i(p), ax^._i(p)q] e Dj, 
7.. : [p,q] eD(X.) - [Yj_,ip),dY:_,ip)q] G D.. 



The proof of Theorem 5.11 does not require the application of the contraction mapping 
principle for nonlinear maps. It uses only Lemma 5.6, Lemma 5.9 and Lemma 5.10 
inductively. 



Proof of Theorem 5.1 1 The main work is the proof of the existence and continuity of X^, 
which we will do first.The assumption of the inductive step is comprised of the following 
two statements. We assume for j: 

1. The map (pi,P2) s D(X^_i) x XjiPi,P2) ^ exists and is continuous. 

Moreover, the pair {x^,X2) — Xj{pi,p2) satisfies 



Xj =£^Pi + QjvIFj_i(xi) 

^2 = ■E]vP2 + Qw■^'^^■fj-l(^1.^2)■ 



(81) 

C82) 



2. The linear map z Qi^Ld^Fj_i{xi,z) maps Dj_^ back into itself and has spectral 
radius less or equal 1/2. 

Both statements of the assumption of the inductive step have been proven for j = 1 in 
Lemma 5.4 and Lemma 5.6 . Let j be smaller than j^^. 



Regularity o1Xj(p) 

Let us first establish that the map 



P = 



Pi 

P2 



■-^ X = 



^;-l(Pl) 

d'Xj_Mp2 
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does not only map continuously into Dj but into for all k < j^^^ — j + 1. 

The argument is the same as in the proof of Lemma 5.8: the map Fj maps continuously 
back into for all k < j^^-j. If x e D| then F/x) e dJ, thus, x = Ef^p+Qf^LFix) e D^+'^ 
for all k < jjnax ~ j (and the dependence on p is continuous because all dependencies are 
continuous) . 



Proof of existence and continuity of d^Xj{p)q 

Let us use the notation p — (Pi,P2) and x — {x^, Xj) = Xj(p). Let Po e D(Xj) be arbitrary. 
We first show that d^Xj{p)q exists for all p in a neighborhood U{pq) with positive distance 
to the boundary of D{Xj). We can choose s > sufficiently small such that p + hqG D{Xj) 
for all h e {—£, e), all q — (^1,^2) with |q| < 1 and all p e f/(Po)- Consider the difference 
quotient 

Xjip + hq)-X^ip) 



1 
h 



Xj_,ip,+hq,)-X^_M 









l(Pl)P2 




Z2 



By assumption of the inductive step, Xj_-^ is continuously differentiable such that the first 
row of this difference quotient has the form 

-1 



Zi(h,pi,qi) = - [Xj_i(pi + hqi)-X^_i(pi)] = d^Xj_^(p^ + hsq^)q^ds 



(83) 



^0 



for h 7^ 0. As established above (Pi,qi) 5^Xj_^(pi)qi e D^ -^ is continuous for all 
k < jnjax ~ J + 1 such that 



Zi(h,Pi,qi)eD;-"'+'^i5f 



2 



(Jmax ~ j + 1 ^2 since j < j^^x), and z{h,Pi,q-^) depends continuously on its arguments, 
also when h = 0. Let us use the abbreviations 

^l(Pl) =^;-l(Pl), 
^2(Pl,P2) = 5%-l(Pl)P2, 

1 

Zi(/i,Pi,qi) = r [^j-i(Pi + ^qi)-^;-i(Pi)] = d^Xj_iipi+hsqi)qids 



JO 



1 



Z2(^,Pl,P2.<?l.<?2) = ^ [^%-l(Pl+^qi)[P2+^<?2] -^^^;-l(Pl)P2] for h / 0. 

With these notations we have Xj_i(pi + hq^) = + hz^ and, for non-zero h, d^Xj_-^(p-^ + 
^<?i)[P2 + ^<?2] = ^2 + hz2. The fixed-point equations (81) and (82) imply a fixed-point 
equation for the difference quotient Z2 for non-zero h: 

Z2 = £^^2 + j^QnL [5 Vj_i(xi + hzi, X2 + hzs) - 5^Fj_i(xi, X2)] 

= £jv<?2 + z(Pi,P2,qi,^) + Qjvf'^^fj-i(^i +^Zi,Z2) where (84) 

^ aiFj_i(xi + fei,X2)-aiFj_i(xi,X2) 
z{Pi,P2,qi,h)^QNL 
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The regularity of Xi, X2 and is: 



Xi e D 



max r^3 

1 — 

1 - ^j-i 



and 



(85) 



We can apply the mean value theorem to the difference quotient appearing in z since x^ 
and z^ are at least in D^_^ and is at least in Dj -^ (see Lemma 5.7, and Lemma A.2 and 
Lemma A.6 in Appendix A) : 



/•I 



z(Pi,P2.qi.^) 



5^Fj_i(xj -\- shzi,X2,Zi,Qi)ds. 



The map (x^, X2,z^,h) ^ 5^Fj_i(x^ +5hz^, X2,Zi, 0) maps Xj, X2, z^ and h continuously into 
the space D^j^l ■' ^ (we see this by applying Lemma 5.7 to d^Fj_i, setting k in Lemma 5.7 

to Jmax - j - !)• Thus, the quantity z(pi,p2, Qi, is in Dj"!^"^ c d]_^ (since j < - !)• 
It depends continuously on p^, (Ji and h in this space, and can be extended to h = 
continuously (such that z{p'^,p2,cii,0) e \ too). 

Hence, (84) is a linear fixed-point problem for Z2 where the inhomogeneity is in Dj"''j ^ 
and depends continuously on {p,q,h). The linear map M(h) : Z2 QjvL5^Fj_i(xi + /iz^)z2 
in front of Z2 on the right-hand side of (84) depends continuously on h as an element of 
LiD^_i',D^_i) (see Lemma 5.10 and note that x^ and z^ are in D^_^). Since the spectral 
radius of the map M(0) (for h = 0) is less or equal than 1/2 by assumption of our inductive 
step, the spectral radius of M(h) is less than unity if we choose h sufficiently small. Thus, 
for all p e D{Xj) and q e E"j and sufficiently small h, Z2 satisfies a contractive linear fixed 
point equation with an inhomogeneity in D|_^ and a contractive linear map that maps into 
D|_j where all coefficients depend continuously on (h,p,q). Consequently, Z2 has a limit in 
D|_j for h ^ that depends continuously on (p, q). For h = the fixed point equation for 



(z^,Z2) simplifies to 



Zi = Ejvqi + Qjv-f'^^-fj-i(^i,Zi) 

^2 = -Ewq2 + Qjv-f' [d^Fj_^ixi,X2,z^,0) + d^Fj_^ix^,Z2)] . 



(86) 



Both equations are linear in q and z = (zi,Z2). Consequently, z(0,p,q), which is by 
definition the directional derivative of Xj in p in direction q, depends linearly on q and 
continuously on p and q. Consequently, 



z(0,p,q) = 



is the Frechet derivative of X, 
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Collection to finish proof of statement 1 of inductive step 

The functions x —Xj{p) and z — d^Xj{p)q satisfy 

X =Ef^p + Qf^LFAx) by inductive assumption (81)-(82), 
z =£jv'3 + QN^^^Pji^'^) by (86) and definition of Fj. 

The variable x — Xj(p) depends continuously on p with respect to the norm of by 
the assumption of the inductive step and the step "Regularity of Xj^p)". The variable 
z — d^X{p)q depends continuously on p and q as shown in the previous step, "Existence 
and continuity ofd^Xip)q'\ Thus (x,z) = iXjip),d^Xip)q) = Xj+^ip,q) e DJ x Dj = D^+i 
depends continuously on (p,q), and satisfies (81)-(82) for j + 1 (which is identical to 
system (87)). This completes the proof of statement 1 of the inductive assumption for 

Spectral radius of map z ^ Qf^Ld^Fj(x,z) 

The map d^Fj maps Dj_,_i continuously into Dj (by Lemma 5.7). Thus, for fixed x the 
linear map z d^Fj{x,z) maps Dj continuously into Dj, and, hence, the map Mj(x) : z 
QffLd^Fj{x,z) maps Dj continuously into Dj, making Mj(x) a compact linear operator. 
Thus, the spectral radius of Mj(x) is determined by its largest eigenvalue (which has 
finite modulus and algebraic multiplicity if it is non-zero). Splitting Mj(x) into its two 
components we get 

QnL a^Fj_i(xi,Zi) 

Q^L [d^Fj_,ix^,X2,z„0) + d^Fj_,ix„Z2)] _ 

If (A, (zi,Z2)) is an eigenpair of Mj{x) then the first row of the definition of Mj{x) implies 
that, either (A,z^) is an eigenpair of z^ — > QjvL5^Fj_i(xi,Zi), or z^ = 0. If (A,Z;^) is 
an eigenpair of z^ Qi^Ld^Fj_i{xi,z-^) then, by inductive assumption, |A| < 1/2. If 
Zi = then the term d^Fj_i{xi,X2,Zi,0) vanishes in the second row, such that (A,Z2) is 
an eigenpair of Z2 Q^Ld^Fj_i{xi,Z2). Thus, by inductive assumption, |A| < 1/2 in this 
case, too. Consequently, the spectral radius of Mj{x) is also less or equal to 1/2, which 
proves statement 2 of the inductive assumption for j + 

Existence of 7j We show inductively that 7j(p) = Fj{Xj{p)). For j — 1 this statement was 
proven in Lemma 5.4. Let j < j^^^^ and assume that Yj — Fj{Xj{p)) for p e D{Xj). Since 

Xjip)^E^p + Q^LFjiXj{p)) 

and Fj maps into D^, Xj is an element of D^. Let q e be arbitrary, and let us denote 
(xi,X2) = iXj{p),d^Xj{p)q) — Xj^i{p,q). The component X2 satisfies 

X2 ^ Ef^q + QNLd^Fjix^,X2) 



Mjix) : 







1 — > 
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such that X2 is in Dj, too. Consequently, 

Yjjp + hq) - _ FjjXjjp + hq)) - F/X/p)) 

h ~ h 

Fjix, + hx2)-Fjix,) FjiXjip + hq))-Fjix,+hx2) ^ ^ 

= Z + Z ■ <;88) 

h h 

Since Fj is continuously differentiable for x^ e D| and deviations hx2 e D| the first quotient 
in the expression (88) converges to d^Fj{xi,X2). Since Fj as a map from Dj into Dj is 
locally Lipschitz continuous the second term in (88) can be bounded by 



F^(Xj(p + hq)) - Fj(Xj(p) + hd'Xj(p)q) 



< 



Xjip + hq)-Xjip) 



-d'Xjip)q 



with some constant K^, for sufficiently small h, which converges to zero for h because 
Xj is differentiable. Consequently, the directional derivative of Yj in p in direction q is 
d^Fj{Xj{p))[dXj{p)q], which is continuous in p and q and linear in q. Therefore, the 
Frechet derivative of 7, exists and 



8-p^^^P\ 



q = d'FjiX^ip),dXj(p)q), 



which implies by definition of Fj and Xj that Yj_^_i — Fj_,_^(Xj+^(p, q)). □ 
We can refine the statement of Theorem 5.11 slightly by noting that Xj : D{Xj) Dj is 
continuous for all j < j^^^ (instead oiXj : D{Xj) Dj). This follows from the continuity 
of Yj = Fj(Xj{p)) as a map into Dj and the relation 

Xjip)^E^p + Q^LYjip). 

Theorem 5.11 completes the proof of the Equivalence Theorem 2.5. The refinement (thatX^ 
maps into Dj) ensures that the image X(p) is in C^'"'^"'"^(T; E"), as claimed in Theorem 2.5 



6. Proof of Hopf Bifurcation Tlieorem 

First, we note that x S{x,co)~^ = x{co~^-) maps C'^(T;]R") into a closed subspace of 
C'^([— T, 0]; M"), if we extend functions x on T to the whole real line by setting x{t) = 
x(tjnod[-7i,7r))- This implies that, if the functional / : C°([-t,0];M") x E M" is EC'' 
smooth then the functional 

(x, pi, co) e C°(T; M") x ^ — /(S(x, co), pt) e M" 

CO 
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is EC'^ smooth, too, such that we can reduce the problem of finding periodic orbits of 
frequency co to the algebraic system (22). The right-hand side Fy in (22) is defined by 

[f (x,co,pi)l (f) = —f(S(A,x,co),^x). 

^ ^ CO 

Let us choose the periodic orbit Xq = (x, co, /i) with x — 0, co — cOq, /i = as the solution 
in the neighborhood of which we construct the equivalent algebraic system. We choose the 
number N of Fourier modes and the size 5 of the neighborhood B°'^(xo) in C°'^(T;M"+2) 
such that the conditions of Theorem 2.5 are satisfied in B^'^{xq). The full algebraic system 
(22) then reads (after multiplication with co and mapping it onto the space rgPj^ from 
^ni2N+i) by applying R-1) 

^PoFyiXyip, CO, ix), CO, ix) + coQoPf^E,^p - Q^P^LFyiXyip, co, pt), co, pt) (89) 

The variables are p e (-^Y^j^id^ -^^as called Py in (22)), and co. We know from 

Theorem 2.5 that 

Yy :(p, CO, /i) e E"(2w+i) >^ ^ X E ^ F{Xy{p, co, /i), co), pi) e C°(T; E"), 
Xy :(p, IX, co) e M"(2jv+i) X K X E ^ Xy{p, co, /x) e C°(T; M") 

are k times differentiable, and note that 

Fy{Xy{0,co,ix),CD),ix) = Q (90) 

for all CO cOq and jit ^ (because Xq = (0, co, pi) is a solution) . The derivative of the 
right-hand side Fy in x — 0, co coq and jx^O with respect to x is A{co,ix)x, defined by 

[Aico,^x)x] (t) = a(M) [xit + co-)], 

where a(/i) is the same linear functional as used in the definition of the characteristic 
matrix K{X,ix) in (25) (the derivatives of F with respect to co and jx are zero due to (90)). 
We observe that A(<D,pt) commutes with Pj and Qj for all j > 0. 

Let us now determine the linearization of Xy{p,co,ii) in (p,co,ix) — (0, co,pi). Due to 
(90) Xy{0, co,ii) is equal to zero for all co ^ coq and ji !^ 0: since is a solution to the 
periodic BVP and PjyO = 0, the zero solution must also be equal to Xy{0, co, jx). Thus, we 
have 



d 

dco 



0= —Xyip,co,^x) 



d 

and 0= — X (p,co,/x) 

p = 



p = 



Moreover, the fixed point equation (77) defining z = [dXy/dp]{p,co,jj,)q, evaluated in 
(p, CO, jx) — (0, CO, jx) reads 

z = £jv<? + QnLA(ix, co) z = £jv<? + Qw^-^Cm, w) Qjv^, (91) 

exploiting that Qj^L = QjvLQjv and QjvA(co,jLi) =A(co,/i)Qjv. In the neighborhood B°'^(xo) 
the spectral radius of Q^LA^fx, co) is less than unity (see Lemma 5.6). Application of Qj^ 
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to (91) gives QjvZ = Q^LA^^, co)Qf^z. Since the spectral radius of QjvLA(/>i, co) is less than 
unity this implies that Qj^z — 0, and, thus 



d 



p = 0. 



q^Ef^q. 



Consequently, the linearization of the algebraic system (89) in (0, co, //) with respect to 
the first variable is 



(92) 



for all <D ^ cOq and /i ^ (also using p for the argument of the linearization in (92)). We 
observe that the linear system (92) decouples into equations for 

yo = Po^nP = -EoP = Po (the average of £^p), 

Yi = QqE^p = p_i sin t + p^ cos t (the first Fourier component of E^^p), 

jj = Qj_iEjp — p_j sin(j t) + Pj cos(jt) (the j-th Fourier component of EjvPj 

2<j< iV), 

where we denote the components of p by pj e E" (j = —N . . .N). This decoupling is 
achieved by pre-multiplication of (92) with Pq and Qj-^Pj for j = 1 . . . AT: 



Po • (92) 
QoPi ■ (92) 
Qj-iPj ■ (92) 



0=A(co,pt)yo = a(M)Po 
= coji -QoIA(to,pi)yi 

0^coyj-QoLAico,^x)yj 



for j = 2...N. 

Inserting the definition of jj into the equations (94) and (95) gives for j >1 



(93) 
(94) 
(95) 



0=co [p_^sin(;t)+P;Cos(jt)] -Qc 



o(/^)[P; sin(;s + jco-) + pj cos(;s + jco-)] ds 



—CO ^p_j sin(jt) + Pj cos(jt)J — : sin(jf )a(/i)[p_j sin(jco-) + Pj cos(jco-)] 
1 

— : cos(jt)a(/x)[— p_,- cos(jto-) + p, sin(j(i)-)]. 
J 

These equations are satisfied if and only if the coefficients in front of sin(jt) and cos(j t) 
are zero. The resulting system of equations reads in complex notation (splitting up again 
into the cases j — 1 and j > 1) 

icouj — a(/i) [u^ exp(i<Ds)] = K(i(D,/i)ui = 0, (96) 
ijcoUj — a(jLt) I^Uj- exp(ijcos)J = K(ijco, jj.) Uj —0 (2 < j < N), (97) 

that is, Uj = p_j + ipj e C" is a solution of (96) (or (97), respectively) if and only if 
— p_jSm{jt) +PjCos{jt) is a solution of (94) (or (95), respectively). 
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The non-resonance assumption of the theorem guarantees that equation (93) is a regular 
Unear system for pg, and that (97) is a regular linear algebraic system for p_j and pj 
ij > 2) at jJL — and CO — cOq (and, hence, for all co and /i near-by). The condition on 
the simplicity of the eigenvalue IcOq oi K ensures that equation (96) (and, thus, (94)) has 
a one-dimensional (in complex notation) subspace of solutions for co — cOq and jJi — 0, 
spanned by the nuUvector of K{ico, 0). Let us denote the adjoint nuUvector of K(icoQ, 0) 
by e C" (again, using complex notation, w"K{icoQ,0) = 0). Since ico^ is simple, the 
relationship 

dK 



wf — (ico,0)vi/0 



holds, which implies that we can choose e C" without loss of generality such that 

r,dK 



w'{—iico,0)v^ = 1- 



With this convention we observe that 



Wi — (ico,0)vi = - — 



/X = 



Cu e C, and 



wf— iC(ico,0)vi = i 
dco 



C (98) 



where Re 7^ by the transversal crossing assumption of the theorem. In complex notation 
any scalar multiple of the nuUvector Vj — v^-\-iv^ is also a nuUvector. Thus, the complex 
scalar factor a + ifS in front of makes up two components of the variable p (in real 
notation) : in short, p solves the linearized algebraic system (92) if and only if all Pj with 
|j| 7^ 1 are zero and p_i sin t + Pi cos t = Re [(a -I- i^)v^ exp(if )] for some a, ^ e E, that 
is. 



"p-l" 


= a 

















(99) 



Let us collect the statements so far and introduce coordinates. We collect all components Pj 
with |j| / 1 and the orthogonal complement in E^" of the space spanned by {bj, into a 
single variable q (of real dimension — n{2N — l)-l-2(n — 1)). Then a set of coordinates 
for p are the variables 



and q e 



We split up the full algebraic system of equations (89) in the same way as we did for the 
linearized problem, by pre-multiplication with Pq and Qj-iPj for j — 1 . . .N: 



^0 

QoPi 



(89) 
(89) 
(89) 



o^PoFiXyip,co,^Ji),co,^Ji) (100) 

= coQoE,p - QoP^LF(Xy(p, co, /i), co, /i) (101) 
= coQj_^EjP - Qj_^PjLF(Xy(p, co, ^jl), co, ^i). (102) 



We split equation (101) further using and its orthogonal complement, the projection 
w| = / — WiW^ /{w^Wi). This gives rise to a splitting into two real equations (w^-(lOl)) 
and 2(n — 1) real equations (w^-(lOl)). Collecting w^-(lOl) and the equations (100) and 



43 



(102) into a subsystem of n{2N — 1) + 2(n — 1) = equations the full algebraic system 
(89) in the coordinates (r,q) has the form 



= 



■M,,(r,q,to,iLi) M,q(r,q,to,pi)' 
Mq,(r,q,to,pi) Mqq(r,q,to,pi) 



By our choice of coordinates the matrices M^^ e E^^^, M^g - ...g^ 
identically zero in r — 0, q — 0, jj, — 0, co — ico^ such that the system matrix has the form 



^2 and M„ 



(103) 



n,x2 





"0 


0" 












"o 


o" 











Mqg(0,0,ico,0) 



(regular) 



(r, q, fi, co) — (0, 0, 0, coq). Thus, we can perform a Lyapunov-Schmidt reduction: we elim- 
inate q by solving the lower equations for q, obtaining a graph q(r, <D,/i) r locally in 
a neighborhood of (r, <d) = (0, 0, 0, cOg). This graph respects rotational invariance: 
q{A^r,co,ii)A^r — Aj[q(r, <D,jLt)r]. Note that the application of to r = (a, ^) corre- 
sponds to the rotation of r by angle s (the same as the multiplication exp(is)(a -I- i^)). The 
Lyapunov-Schmidt reduction of (103), replacing q by the graph q{r, co,ijl) r, then reads 



= M^j.{r, q(r, co, /i) r, co, jj,) r —: M^{r, co, /i) r. 



(104) 



where M,. is still rotationally symmetric in r: Mr(A5r, <D,jLt) A^r = A^M^r, co,ijl) r. Equa- 
tion (98) in real notation implies that 



—^(0, cOo, 0) = ^(0, 0, cOo, 0) = 
o CO a CO 

dM, dM„ 

—^(0, cOo, 0) = ^(0, 0, cOo, 0) = 



-1 

1 

Rec^ 
Imc 



Imc^ 
Rec 



Equation (104) is a system of two equations with four unknowns (r = (a, (3), co and /i). 
We now fix one of the unknowns setting 

a = 

such that we can expect one-parametric families of solutions (^, co,pi). Introducing Mp 
as the second column of and dropping the dependence on a (which is zero), the first 
derivative of MJjd, co, jj.) in (0, cOq, 0) with respect to the pair co and jJL is: 



dco d jJL 



(0,cOo,0) 



-1 




Imc^ 
Rec 



which is regular (as established in (98), since Rec^ 7^ due to the assumption that the 
eigenvalue crosses the imaginary axis transversally) . Note that M« itself is a projection of 
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the first derivative of the original right-hand side of the full algebraic system (89). Thus, 
Mp is — 1 times continuously differentiable, and we end up with a system of two equations 
for three scalar variables (13, co, ji): 

Hence, either 13 — 0, which corresponds to the trivial solution or (after division by (3) 

0^Mp{f3,co,lx), (105) 

where M^(0, (Dq, 0) = (0,0) and the derivative with respect to the pair (<x),pt) is regular 
in (0, cOojO)- Thus, we can apply the Implicit Function Theorem to (105) to obtain a 
unique graph (to(^),jLi(^)) solving (105). The graph satisfies ((d(0), /i(0)) = (tOg, 0), 
and, thus, branches off from the trivial solution (which has 13 = and co and jj, arbi- 
trary) . The rotational symmetry of M,. implies reflection sjmimetry of Mp in (3 such that 
Mp{—I3, CO, ijl) — Mp{l3, CO, jj.) for all 13, co and jJL. Hence, the solution graph is reflection sym- 
metric, too: co(— ^) = co(^) and iLi(— ^) = jLt(^). Thus, for small ^ there is a unique non- 
trivial solution of the full algebraic system of the form r = {0, (3), q = q(r, co(^), jLt(^)) r. 
As Equation (99) shows, we can extract the coordinates a (which is zero) and (3 from the 
full solution X e C'^(T; E") by applying the projections 



1 

71 
1 

71 



" 71 



cos(t)vJx(f ) - sm(t)vjx(t) dt = 



71 



71 



71 
" 7T 



sin(f )v^x(f ) -I- cos(f )v^x(t) dt = — 

'' ' 71 



Re exp(it)] ^x(t)dt = a, 
Im [viexp(it)]^x(t)df = 



which determines the First Fourier coefficients of x as claimed in (27) in Theorem 3.2. 
(Recall that the vector = v^-l-V; was scaled to have unit length and that the decomposition 
was orthogonal.) □ 
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A. Basic differentiability properties of tlie riglit-liand side 

Let J be a compact interval or T. Let (D, || • H^) be a Banach space of the form 

D = C''i(J;E'"i) X ... X C^t(J;R"'t) 

where i >1, the integers kj are non-negative and the integers mj are positive. We use the 
natural maximum norm on the product D: 

\\x\\D = Kxi,---,Xti)\\D= max ||x ||k, 

Mi,...,t] ' 

and use the notation 

D'< = C^^+\J;W^i) X ... X C'^^+'^(J;M'"0, \\x\\Dk = max \\x^^\, 

0<j<k 

D°'^ = {x eD : L(x) < oo}, with the norm = max{||x||£„L(x)}, 

where x^^^ is the component- wise jth derivative and the Lipschitz constant L(x) is defined 
as 

\xf'\t)-xf'\s)\ 

L(x) = sup max ; ; , 

t#s J=W \t-s\ 



47 



where t and s in the index of sup are taken from J, if J is a compact interval, and from E 
if J = T. Balls that are closed and bounded in D°'^ are complete with respect to the norm 
of D. 



A.1 . Basic properties of / 

This section proves three properties that EC^ smooth functionals / have: first that the 
derivative limit (10) exists also for Lipschitz continuous deviations, second a weaker form 
of the mean value theorem, and third, local EC Lipschitz continuity. 

Lemma A.1 (Extension of derivative to deviations in D^'^) 

Let / : D "-^ M" be EC^ smooth in the sense of Definition 2.1. Then the limit 
required to exist in Definition2.1 exists also in the || • H^^-norm: for all x e 

r l/(^+y)-/W-gV(x,y)| 

lim 7—r. = 0. (106) 



Note that in (106) the norm in which y goes to zero is || • ||q ^ instead of || • \\d i- 



Proof This is a simple continuity argument. Let £ > be arbitrary. We pick 5 > such 
that 

l/(x +y) -/ W - 5 V(x,y)| < sWyWd.i (107) 

for all y satisfying ||y ||d i < 5. Let y e D°'^ be such that ||j||o,l < ^- can choose a 
y e that satisfies 

||y||o,i<min{5,2||y||o_a (108) 

\fix + y)-fix + y)\<^\\y\\D,L (109) 
\d'f(x,y-y)\<e\\y\\^^,. (110) 

Condition (108) can be achieved because is a dense subspace in D°'^, and for every 
element y of the || • 11^,^-norm is not larger than the || • H^j^^-norm: Hyllo i < ||y|lD_L- 
(109) follows from the continuity of / and the density of D°'^ in D\ and (110) follows 
from the continuity of d^f as a map on x D, and the density of D°'^ in D^. Combining 
estimate (107) with (108)-(110) we obtain 

l/(x + y) -/(x) - a V(x,y)| < 4^||ylUi. 

□ 

Lemma A.2 (Existence of mean value) 

There exists a continuous function 

I a-.D^xD^xD^R" 
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I 



which is hnear in its third argument and satisfies for all x,y gD^ 
Moreover, a(x,0,y) = d^f{x,y) for all x G and y 



(Ill) 



Proof The argument for the existence of a mean value follows exactly the proof of the 
general mean value theorem [12] : the candidate for d{u, v, w) is 



a(u, v,w) = 



d^f{u +sv,w)ds. 



(112) 



Since d is assumed to be continuous in its arguments the integral is well defined and 
continuous in its arguments ugD^,vgD^,wgD. All one has to show is that the a 
defined in (112) satisfies (111): let x,y G and £ > be arbitrary, and choose m such 
that uniformly for all s e [0, 1] 



1 k \ 

d'fix+sy,y)ds--yd'f{x + -y,y 

f {^+'y + ^) -fix+sy)-d'f (^x+sy,^^ 



< 



m 



Then it follows that 



/(x + y)-/(x) 



d'^f{x+sy,y)ds 



< 2s. 



Since e > was arbitrary the left-hand side must be zero. □ 

Lemma A.3 (Local EC Lipschitz continuity) 

For all X e D°'^ there exists a neighborhood U(x) c and a constant > 
such that for all y^ and 72 ^ U{x) the following Lipschitz estimate holds: 

l/(yi)-/(y2)l<^xllyi-y2llD- 



Note that the upper bound depends only on the || • 11^,-norm, not on the || • |Id,l 
would be a weaker statement. 



norm, which 



Proof We prove the Lipschitz continuity first for y^ and 72 from a sufficiently small 
neighborhood [/(x) n c of x e 

Let X be an element of D^. Since the mean value a is continuous in (x, 0, 0), and 
a(x, 0, 0) = 0, we have a 5 > such that for all u,v gD^ and w G D satisfying j < 5, 
\W\\d,i < 5 ^i^d ll^llo < 5 

|a(x -l-u,v,w)| < e. 
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This implies that \d(x+u, v, w)| < [e/5]||w||o for u and v with max{||u||o j, HvH^, 1} < 5 and 
w G D (since a is hnear in its third argument) . Thus, ||a(x + u, v, •)|lD<^/5fora(x + u, v, •) 
as an element of L(D;D) in the operator norm corresponding to D. Consequently, if 
llyi - ^IId.i < 5/2 and - ^IId,i < 5/2 



l/(yi)-/(y2)l = 



5(72, Ji - 72,71 -72) ds 



<dlyi 



such that we can choose K^ = e/5. The extension of the statement to D°'^ follows from 
the continuity of / in D: U{xq) n is dense in U{xq) c D°'^ using the || • H^^-norm. Pick 
two sequences y„ and z„ in U{xq) n that converge to y and z in U{xq] in the Lipschitz 
norm. Then fiy„) — > /(y) and /(z„) fiz) since / is continuous in D. Moreover, 
lljn -zJId Wy-zWo for n ^ 00. Since 



\fiyn)-fiZn)\<Kj\yn-zJo 

for all n the inequality also holds for the limit for n ^ 00. 



(113) 

□ 



A.2. Basic properties of F 

In this section we restrict ourselves to the periodic case: J = T. Let F : D ^ C°(T; M") be 
defined as F(x)(f) =/(AfX). 

Lemma A.4 (Continuity of F) 

Let / : D E" be continuous. Then F : D ^ C°(T; E") is also continuous. 



Proof This is a simple consequence of the continuity of/, the continuity of(f,x)-^AfX 
with respect to both arguments (t and x) in the || • ||o-norm, and the compactness of T. Let 
£ > and X e D be arbitrary. We want to prove continuity of F in x. So, we have to find a 
5 > such that 

|/(A,x + h)-/(A,x)| < e for alls eT and he D, satisfying ||h||o < 5. (114) 

(Since || AjhH^ = HhH^, we can replace A^h by h.) The continuity of / implies that for every 
r > and every t e T we find a 5^{t, r) such that 

|/(AtX + h)-/(AfX)| < r whenever ||h||o < 5^(t,r). (115) 

For every t e T there exists an open neighborhood [/(t) c T such that 

||A,x - AfXilo < 5^(t, e/2)/2 for all s e [/(t), 

because the function t e T A^x is continuous in t. These neighborhoods U{t) are an 
open cover of the compact set T, so there exist finitely many t^, . . . , e T such that the 
union of the neighborhoods U{tj) contains all points s e T. We choose 

5= min 5,(tj,£/2)/2, 

j=l,...,m 
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which is a positive quantity. Let s e T be arbitrary and let h e D satisfy ||h||£, < 5. We have 
to check the inequahty (114). The point s is in one of the neighborhoods U{tj), say without 
loss of generality, s e [/(ti)- Thus, HA^x — Aj^xH^, < 5^(ti, e/2)/2, and, consequently, 
llAjX - A(.^x + h||o < 5^(ti,e/2) (because also HhH^, <5 < 5^{t-^, le/l)/!). Therefore, we 
can split up the difference \f{A^x + h) — /(A^x)!: 



|/(A,x + h)-/(A,x)| 



< 



[/ (a,^x + (A,x - A,^x + h)) - /(A,^x)] 
[/ (a,^x + (A,x - A,^x)) -/(A,^x)] 



+ 

< e/l + e/l = e 

Note that the deviations from A^^x in the arguments of / in both terms of the sum are less 
than or equal to 5^{ti,el2) such that we can apply (115) for t — t-^, r — e/2. □ 
The following lemma lists properties that F has if / satisfies local EC Lipschitz continuity 
in the sense of Definition 2.2. That is, we do not assume that / is EC^ smooth in the sense 
of Definition 2.1 for Lemma A. 5. Since Lemma A. 3 was proved using only the assumption 
of EC^ smoothness of /, local EC Lipschitz continuity is a weaker condition. 

Lemma A.5 {EC Lipschitz continuity of F) 

Assume that / : D — > E" is locally EC Lipschitz continuous in the sense of Defini- 
tion 2.2. Then F has the following properties: 

1. for all X e D°'^ there exists a neighborhood f/(x) c and a constant 
K^>Q such that for all and 72 ^ U{x) 

\\F{yi)-F{y2)\\o<K\\yi-y2\\D- 



2. F maps elements of D°'^ into C°'^(T;]R"). Moreover, for every x e D°'\ any 
bounded neighborhood f/(x) c for which the Lipschitz constant K^. 
exists has a bounded image under F : there exists a bound R> such that 
ll^^(y)llo,i < R for all y e L/(x) (R depends on [/(x)). 



Proof Statement 1 is a consequence of the local EC Lipschitz continuity of / and the 
compactness of T (which allows one to choose a uniform Lipschitz bound for all t e T) . 

Concerning statement 2: let x e D°'^ be arbitrary, and let the neighborhood U{x) be 
bounded (say, ||y — xH^, ^ < 5) such that F has a Lipschitz constant in U{x). Then we 
have for all y,z e U{x) and f,s e T the estimate 

|/(A,y) -/(A,z)| < Kj\A,y - A.zWd = ^xllA,_,y -zW^. 

Initially setting z = x ands = t we get abound on ||F(y)||o: ||F(y)||o < ||F(x)||o+i<^x^ —'■ ^0 
for all y e L/(x). It remains to be shown that the Lipschitz constant of F(y) is bounded for 
jeL/(x): 

\FiyXt)-FiyXs)\ = \fiA,y)-fiA,y)\<Kj\A,y-A,y\\,,<Kj\y\\j,Jt-s\. 
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Since \\y — x\\q i < 5 for y e U(x), choosing 

R^max{Ro,K,(\\x\\o^, + 5)} 

ensures that ||F(y)||o ^ <R. □ 
Define the maps 

a^F(u,v)(t) = d^fiA,u,A,v) for ugD^,vgD, 

A{u, V, w)(f ) = a(AtU, AfV, A^w) for u e D\ v e D\ w e D. 

The following Lemma A. 6, and Lemma A. 7 assume that / is £C^ smooth in D in the sense 
of Definition 2.1. 

Lemma A.6 (Differentiability of F) 

Let / : D — > E" be EC^ smooth. Then F, d^F and A have the following properties. 

1. The map (u, v) 5^F(u, v) is continuous in both arguments (and linear in 
its second argument) as a map from x D into C°(T; E"). It satisfies for all 

XGD^ 

\\Fix + y)-F{x)-d'F{x,y)\\, 
hm — =0. (116) 

y^D"^' WyWd.l 

2. The map A(u,v,w) is continuous in all three arguments (and linear in its 
third argument) as a map from x x D into C°(T; M"). It satisfies for all 
x,yGD^ 

F(x + y)-F(x)^A^(x,y,y). 
Moreover, A(x,0,y) = 5^F(x,y) for all x G and y e D. 

Note that in the limit (116) we allow for deviations y e D^'^. 

Proof The continuity of d^F follows from the continuity of d^f by applying Lemma A.4 
to d^f : X D >— > E" instead of/. The linearity of d^F in its second argument follows 
from the linearity of d in its second argument. 

The limit (116) also follows from the corresponding limit (106): let x G and e > 
be arbitrary. For every fixed t there exists a 5(t) > such that 

|/(A,x + A,y) -/(A,x) - a V(A,x, A,y)| 

<s (117) 



\D,L 

for all y with ||y||oi < 5(f). As / and d^f are continuous in their arguments x G and 
y G D°'^, the inequality also holds for all s in a sufficiently small open neighborhood of t, 
U(t). The set of neighborhoods U{t) for all t e T are a cover of the compact set T by open 
sets. Choosing a finite subcover from this cover, and labeling the times f i, . . . , t^, we can 
choose 

5 — min 5(t^) 

k=l,...,m 
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such that (117) holds for all uniformly t e T. This proves statement 1 of the lemma. 

Concerning statement 2: for the continuity of A we invoke again Lemma A.4, this 
time for a onD^ x x D. The linearity of A in its third argument follows from the 
linearity of a in its third argument. The relations F(x + y)(t) — F(x)(t) = A{x,y,y){t) 
andA(x,0, y)(f) = d^F{x,y){t) in every t e T follow from the corresponding relations for 
/ and d, as stated in Lemma A.2. □ 

Lemma A.7 (Differentiability of images of F) 

Let / : D E" be £C^ smooth and let > be some integer. We assume that 
F -.D^ C''(T;]R") and d^F : x D ^ CHT;E") are continuous maps. Then F 
maps elements of into C'^^^(T; E"), and F is continuous as a map from to 
I C'^+HT;M"). 



Proof Let x be in D\ that is, x' gD. If d^F : x D ^ C''(T;M") is continuous then 
A : X X D ^ C^T;^), which is given by A(u,v,w) = d^F^u + sv,w)ds, is 



continuous, too. Using statement 2 of Lemma A. 6 we have 



F(A,x)-F(x) 



A X, Aj^x — X, 



(118) 



On the right side || A^jX — x||o ^ converges to for h — > 0. Also, 



for h ^ 0, 



because x e D\ Since A is continuous in its arguments the right side converges to 
A(x,0,xO = 2^F(x,xO e C''(T;E") for h ^ 0. Thus, the limit of the left-hand side in 
(118) for h ^ exists, too, such that F(x) e C''+^(T; E") and the time derivative (F(x))' 
equals d^F<ix,x'). Since (v,w) G x D ^ 5^F(v,w) e CHT;E") is continuous in (u,v), 
the time derivative of F(x), (F(x))' = d^F{x, x') is also continuous in x if we use the norm 
II • ll^, 1 for the argument and || • for the image. □ 
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